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NOMENCLATURE 


A  dimensional  constant;  activation  energy 

a  wire  radius;  dimensional  constant 

B  dimensional  constant 

b  dimensional  constant 

C  vapor  concentration 

Cp,  c  heat  capacity 

Cj  fractional  increase  in  snow  density  per  meter 

water  equivalent  of  load  per  hour  at  ps  =  0 
and  0  =  0°C 

C2  dimensional  constant  of  compaction  par¬ 
ameter  of  snow  (m3/mg) 

C3  dimensional  constant  (fractional  settling 

rate  at  0°C) 

C4  dimensional  constant 

D  diffusion  coefficient 

F  load  on  the  object 

F  force  per  unit  length 

G  air  mass  flow  rate 

H  enthalpy 

h  thickness  of  a  liquid  layer 

L  latent  heat 

M  molecular  weight;  molecular  weight  of  air 

m  net  sublimation  per  unit  volume;  mass;  mass 
flux 

P  hydrostatic  pressure;  vapor  pressure 

APq  maximum  excess  pressure 

p  pressure;  precipitate 

R  gas  law  constant 

5  fractional  salt  content;  entropy  per  unit  mass 

t  time 

T  temperature  in  kelvins 

U  internal  energy 

V  volume  of  an  object,  ice;  fractional  volume 

content 

v  penetration  velocity;  volume  per  unit  mass 

weight  of  snow  above  a  given  layer  expressed 
in  terms  of  water  equivalent  in  meters 

w  mass  of  unfrozen  water 

a  thermal  diffusivity;  proportional  constant 

fa  ratio  of  vapor  concentration  to  temperature 

7  coefficient  of  thermal  expansion 


6  thickness  of  coating;  water  layer 
r)  viscosity  coefficient  of  snow 

viscosity  coefficient  of  snow  when  snow  density 
is  extrapolated  to  zero 

0  temperature  in  degrees  Celsius 

A  thermal  conductivity 

n  3.1 41 6;  total  pressure 

p  density;  thermal  resistivity 

a  salinity  (grams  of  salt  per  gram  of  sea  ice) 
r  time  increment 

0  correlation  coefficient 

cj  compressibility 

Subscripts 
a  air;  adiabatic 

b  brine 

bi  bubbly  ice 

c  cubic 

e  effective;  experimental 

f  fusion 

h  precipitated  hydrate 

i  ice 

ia  ice  containing  air  bubbles 

£  liquid;  linear;  linear  thermal  boundary 

m  melting;  material;  mean 

n  nth  layer 

o  through  air 

p  pressure 

s  snow;  solid;  sinusoidal  boundary;  sublimation; 

saturation;  and  distance  increment 

sa  snow  containing  air 

se  snow  (effective) 

si  sea  ice 

t  isothermal,  theoretical 

v  volume;  vapor 

w  water 
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REVIEW  OF  THERMAL  PROPERTIES  OF 
SNOW,  ICE  AND  SEA  ICE 


Yin-Chao  Yen 


INTRODUCTION 

This  review  was  undertaken  in  an  attempt  to  sum¬ 
marize  and  analyze  as  completely  as  possible  the  re¬ 
ported  data  on  the  various  thermal  properties  of 
snow,  ice  and  sea  ice  and  to  provide  readily  available 
information  for  practical  use  in  the  field  of  snow  and 
ice  research. 

In  the  second  section,  research  on  the  density, 
linear  and  cubic  expansion  coefficients  and  compres¬ 
sibility  of  ice  is  summarized.  Slight  variations  in 
density  due  to  various  defects  such  as  contamination 
(by  solids  as  well  as  by  air),  aging,  and  crystallographic 
form  are  discussed.  In  the  temperature  range  of  prac¬ 
tical  interest  in  cold  regions,  both  the  linear  and  cubic 
coefficients  of  expansion  can  be  satisfactorily  ex¬ 
pressed  as  linear  functions  of  temperature.  The  re¬ 
ported  '  ork  m  isothermal  and  adiabatic  compressi- 
bilit)  a tii  limited  and  the  variation  of  results  of 
different  m  stigators  is  pronounced. 

In  the  tl  d  section,  snow  density  changes  due  to 
various  metamorphism  processes  are  discussed.  These 
phenomena  are  important  because  of  their  effect  on 
the  physical  and  mechanical  properties  of  snow  and 
on  the  physical  processes  occurring  within  a  snow 
mass. 

The  process  of  regelation  is  described  in  the  fourth 
section.  A  great  number  of  theoretical  analyses  and 
experimental  works  are  thoroughly  reviewed  and 
compared.  Discrepancies  between  the  proposed  the¬ 
ories  and  the  experimental  results  are  great.  The  re¬ 
lationships  between  the  penetration  velocity  of  the 
object  passing  through  ice  and  the  wire  material, 
wire  size,  stress  level  and  ice  purity  are  also  discussed. 
Knowledge  of  the  physics  of  regelation  may  provide 
an  insight  into  the  processes  of  sintering  and  the  de¬ 
velopment  of  intergranular  bonds  in  snow. 

In  the  fifth  section,  the  heat  capacity  and  ther¬ 
mal  conductivity  of  snow  and  fresh-water  ice  over  a 


great  range  of  temperatures  are  reviewed,  and  expres¬ 
sions  developed  by  regression  analysis  are  presented. 
These  properties  control  the  rate  of  propagation  of 
thermal  waves  through  the  snow  and  ice  mass  and 
indicate  both  the  relative  potential  as  a  heat  storage 
medium  and  the  rate  of  heat  dissipation.  The  effects 
of  water  vapor  diffusion  (under  either  natural  or 
forced  convection)  on  the  temperature  profile  and 
mass  redistribution  within  a  snow  layer  are  discussed, 
providing  a  basis  for  interpreting  the  field  data. 

Finally,  a  general  discussion  of  the  thermal  prop¬ 
erties  of  sea  ice  is  given.  Since  sea  ice  is  a  much  more 
complicated  material  to  deal  with  than  pure  ice,  and 
its  composition  (the  relative  proportion  of  its  con¬ 
stituents)  is  strongly  temperature-  and  time-dependent, 
few  experimental  studies  have  been  reported.  It  is  be¬ 
lieved  that  predicted  thermal  properties  of  sea  ice, 
based  on  some  simplified  sea  ice  structural  models, 
give  a  good  approximation  of  the  real  value.  In  this 
review,  expressions  relating  specific  heat,  heat  of  fu¬ 
sion,  density,  thermal  conductivity  and  air  bubble 
content  for  some  specific  temperature  ranges  are  given. 
Thermal  conductivity  models  and  methods  of  deter¬ 
mining  thermal  diffusivity  are  also  briefly  described. 

DENSITY,  THERMAL  EXPANSION 
AND  COMPRESSIBILITY  OF  ICE 

Density 

According  to  a  review  work  by  Dorsey  (1940),  the 
bulk  density  of  ice  at  0°C  and  atmospheric  pressure 
varies  from  0.916  to  0.918  Mg/m3.  Barnes  (1901) 
and  Dantl  and  Gregora  (1968)  indicated  that  ice 
density  pt  decreases  slightly  with  age.  Nichols  (1899) 
reported  densities  at  0°C  of  0.91795,  0.91632  and 
0.91603  Mg/m3  for  freshly  formed  natural  ice,  one- 
year-old  natural  ice  and  artificial  ice  frozen  at  low 
temperatures,  respectively.  However,  Barnes  found 
much  smaller  variation  in  the  ice  samples  taken  from 
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Figure  1 .  Density  of  ice  as  a  function  of  temperature.  A-  deduced  by  Lonsdale 
(1958)  from  the  average  of  several  diffraction  measurements  made  prior  to  1958. 
o -deduced  by  Eisenberg  and  Kauzmann  (1969)  from  X-ray  diffraction  measure¬ 
ments  by  LaPiaca  and  Post  (I960).  •— deduced  by  Hobbs  ( 1 974)  from  X-ray 
diffraction  measurements  made  by  Brill  and  Tippe  (1967). 


the  St.  Lawrence  River;  densities  at  0°C  of  new,  one- 
year-old  and  two-year-old  ice  were  0.91662,  0.91648 
and  0.91637  Mg/m3,  respectively. 

There  are  many  other  factors  affecting  the  varia¬ 
tion  of  the  bulk  density  of  ice  pj.  These  include  the 
number  and  nature  of  the  cracks,  the  degree  of  air 
entrainment,  ice  purity,  dislocation,  and  stacking 
fault  vacancy. 

The  true  value  of  pv  reported  to  be  0.9167 ± 
0.00005  Mg/m3,  was  determined  by  Ginnings  and 
Corruccini  (1947)  using  a  Bunsen  ice  calorimeter;  pj 
may  also  be  deduced  from  measurements  of  the  unit¬ 
cell  parameters.  Figure  1  shows  pj  versus  1 1T  from 
data  obtained  by  Lonsdale  (1958),  Eisenberg  and 
Kauzmann  (1969)  and  Hobbs  (1974).  The  value  of 
Pi  at  0°C  given  by  Lonsdale  (0.9164  Mg/m3)  is  in 
good  agreement  with  that  given  by  Ginnings  and 
Corruccini.  However,  for  lower  temperatures,  pj 
values  given  by  LaPiaca  and  Post  (1960)  and  Brill 
and  Tippe  (1967)  are  probably  more  reliable. 

Thermal  expansion 

The  coefficient  of  linear  expansion  is  a  measure 
of  the  fractional  change  in  length  per  unit  change  in 
temperature.  Butkovich  (1957)  reported  that  the 
orientation  of  the  C-axis,  the  type  of  ice  (whether 
single  or  polycrystalline),  and  the  grain  size  do  not 
appreciably  affect  the  values  of  the  coefficient  of 
linear  expansion.  Ice  can  be  considered  as  an  iso¬ 
tropic  material  with  respect  to  thermal  expansion  in 
the  temperature  range  0°  to  -30°C.  Figure  2  shows 
some  of  the  most  reliable  measurements  of  7g  for 
bulk  ice;  it  can  be  seen  that  7g  increases  with  increas¬ 
ing  temperature.  According  to  Jakob  and  Erk  (1928) 


and  Dantl  (1962),  7g  is  negative  at  about  70  K  and 
lower.  Hamblin  (Powell  1958)  reported  that  7g  values 
measured  in  a  direction  parallel  to  the  C-axis  are  about 
1.8%  and  10%  greater  than  those  in  a  direction  perpen¬ 
dicular  to  the  C-axis  at  273.1  and  73.1  K,  respectively. 


Figure  2.  Coefficient  of  (inear  expansion  of  ice  at 
atmospheric  pressure.  Points  represent  the  data  of 
Jakob  and  Erk  (1928),  Powell  (1958),  Butkovich 
( 1 959)  and  Dantl  ( 1 962).  Samples  consisted  of  poly¬ 
crystals,  single  crystals  and  single  crystals  parallel  and 
perpendicular  to  the  C-axis.  (The  data  of  Dantl  are 
taken  from  his  graph  and  therefore  may  introduce  a 
slight  error.) 
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Figure  3,  Coefficient  of  linear  expansion  of  ice  de¬ 
duced  from  measurements  of  the  temperature  de¬ 
pendence  of  the  unit-cell  parameters . 


Dantl  found  that  D20  ice  has  a  slightly  higher  value 
for  throughout  the  temperature  range  studied.  His 
7C  values  at  240,  260  and  270  K  are  believed  to  be  in 
error.  For  temperatures  T  greater  than  80  K,  the  dak 
can  be  represented  by 

■ygx  10*  =  -1 1.7582  +  0.2424  T  (1) 

with  a  relatively  higher  correlation  coefficient  <j>  of 
0,9736.  Therefore,  for  engineering  applications,  it 
is  quite  adequate  to  compute  7C  from  this  expression, 
which  is  applicable  in  the  temperature  range  from  80 
to  273.1  K. 

!n  Figure  2  the  data  points  of  Powell  were  calcu¬ 
lated  from  two  expressions: 

7fi  x  1 06  =  56.5  +  0.250  0  (2) 

for  measurements  in  a  direction  parallel  to  the  C-axis 
and 

7fi  x  1 06  =  55.5  +  0.248  6  (3) 

for  measurements  in  a  direction  perpendicular  to  the 
C-axis,  where  B  is  temperature  in  degrees  Celsius. 

Values  for  7ft  can  also  be  deduced  from  measure¬ 
ments  of  the  temperature  dependence  of  the  unit-cell 
parameters  of  ice.  Values  obtained  this  way  may  be 
more  significant  that  those  obtained  by  direct  diato- 


Figure  4.  Coefficient  of  cubic  expansion  of  ice  as  a 
function  of  temperature  at  atmospheric  pressure, 
o- deduced  by  Lead  belter  (1965)  from  linear  expan¬ 
sion  coefficients  given  by  Powell  (1958)  and  Dantl 
( 1962).  •— deduced  by  Eisenberg  and  Kauzmann 
(1969)  from  measurements  of  the  unit-cell  parameters 
made  by  LaPlaca  and  Post  ( 1 960).  A  -  deduced  by 
Hobbs  (1974)  from  measurements  of  the  unit-cell 
parameters  made  by  Brill  and  Tfppe  (1967). 

metric  measurements  because  they  depend  only  on 
changes  in  the  dimensions  of  the  lattice  and  not  on  the 
texture  of  the  ice.  However,  X-ray  measurements  may 
involve  larger  errors  than  bulk  measurements  at  low 
temperatures.  Figure  3  shows  7C  values  deduced  from 
the  temperature  dependence  of  the  unit-cell  parameters 
from  more  recent  measurements  made  by  LaPlaca  and 
Post  (1960)  and  Brill  and  Tippe  (1967).  The  magni¬ 
tude  and  general  trend  seem  to  be  in  fair  agreement 
with  the  measurements  for  bulk  ice  shown  in  Figure 
2.  For  temperatures  from  80  to  273  K,  7C  can  be 
fairly  represented  by  7cx  106  =  -15.48  +  0.28  T.  There 
seems  to  be  no  consistent  difference  between  7C  values 
parallel  and  perpendicular  to  the  C-axis  in  either  set  of 
measurements.  The  irregular  pattern  of  the  data  by 
LaPlaca  and  Post  must  be  due  to  a  lack  of  sample 
purity  or  to  experimental  conditions.  For  practical 
purposes,  eq  1  can  be  used  to  calculate  7C  values  for 
temperatures  ranging  from  80  to  273  K. 

Figure  4  shows  the  coefficient  of  volumetric  ex¬ 
pansion  7C  of  ice  as  a  function  of  temperature  at  at¬ 
mospheric  pressure.  The  data  were  taken  from  Hobbs 
(1974)  [omitting  two  7C  values  at  mean  temperatures 
of  163.1  and  183.1  K  from  measurements  of  the  unit¬ 
cell  parameters  made  by  Brill  and  Tippe  (1967)) .  A 
linear  regression  analysis  results  in  an  expression  of 
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7cx  106  =  -24.86+0.67  T  with  a  fairly  high  correla¬ 
tion  coefficient  of  0.94.  Though  there  are  some  dis¬ 
crepancies  among  the  three  sets  of  data,  the  general 
trend  of  yc  variation  with  temperature  is  evident:  It 
increases  with  temperature  and  it  becomes  negative 
when  t  is  lower  than  about  50  K. 

Compressibility 

The  compressibility  of  a  substance  is  defined  as 
the  change  in  its  volume  per  unit  change  in  hydro¬ 
static  pressure.  If  the  change  takes  place  at  constant 
temperature,  it  is  called  isothermal  compressibility 
c*;t.  If  the  change  takes  place  without  energy  exchange 
with  the  surroundings,  it  is  called  the  adiabatic  com¬ 
pressibility  u;a. 

As  usual,  the  limited  data  available  did  not  agree. 
Bridgman  (1912)  reported  an  oot  of  37 x  10“6  /bar 
at  273.1  K  and  one  atmosphere;  this  value  is  about 
three  times  higher  than  the  value  12x  1CT6  at  70°C 
and  300  atmospheres  reported  by  Richard  and  Spey- 
ers  (1914).  Bridgman  later  revised  his  values  as 
shown  in  Figure  5,  which  indicates  the  effect  of  T 
on  ojv  With  the  exception  of  u>t  at  273  K,  the  five 
points  given  by  him  lie  close  to  a  line  on  a  semi-log 
plot.  The  values  for  cut  can  be  given  by  cutx  1 06  = 
14.20  exp  (0.0018  T). 

Values  for  cua  obtained  by  L eadbetter  (1965) 
based  on  the  elastic  constant  measurements  of  ice  by 
Bass  et  al.  (1957)  and  Zarembovitch  and  Kahane 
(1964)  are  also  shown  jn  Figure  5  and  can  be  expressed 
as  cuax  106  =  10,55  exp  (0.0007  T).  Leadbetter  indi¬ 
cated  that  the  uncertainty  in  coa  values  for  tempera¬ 
tures  below  -30°C  is  probably  less  than  10%.  For  tem¬ 
peratures  above  -30°C,  u>a  values  are  accurate  to  about 
5%.  Dantl  (1969)  also  presented  cja  values  based  on 
measurements  of  the  elastic  constant  and  provided  the 
following  expression: 
coax106  =  11.94 

(1+1 .653x  1 0"3  0  +  3.1 2x  1 0'6  d 2 )  (4) 
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Figure  5.  Isothermal  and  adiabatic  compressibility  of 
ice.  A— data  of  Bridgman  ( Richard  and  Speyers  1914). 
o- deduced  by  Leadbetter  (1965)  from  elastic  com 
stants  measured  by  Bass  et  al.  (1957)  and  Zarembo¬ 
vitch  and  Kahane  ( 1 964).  • —data  from  Dantl  ( 1 969). 

DENSITY  CHANGES  IN  SNOW 
Compaction 

The  compaction  of  snow  layers  can  be  seen  graph¬ 
ically  in  seasonal  snow  layer  depth  profiJes  obtained 
during  several  investigations  (Bader  et  al.  1939).  Ko- 
jima(1967)  and  Yosida  (1963)  presented  quantitative 
expressions  for  density  change  due  to  compaction. 
Based  on  many  observations  of  the  change  in  depth 
of  various  layers  (with  no  change  in  water  equivalent), 
the  following  relation  between  snow  density  and  over¬ 
burden  weight  of  snow  was  established  : 


where  cua  is  in  bar'1  and  0  is  in  degrees  Celsius.  A 
few  points  from  this  expression  are  also  shown  in  the 
figure.  Note  that  Dantl’s  o?a  values  are  much  lower 
than  Leadbetter’s  values  and  that  they  decrease 
sharply  as  temperatures  decrease.  According  to 
Dantl’s  data,  wa  becomes  negative  around  120  K. 

On  the  other  hand,  Leadbetter  found  that  the  varia¬ 
tion  of  cua  with  T  was  very  slight.  This  discrepancy 
must  be  due  to  a  variety  of  factors,  for  example,  the 
conditions  of  the  experiment  and  the  origin,  age  and 
purity  of  the  sample. 


Ps  df  i? 
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where  ps  is  the  snow  density  (Mg/m3),  t  is  the  time 
(hr),  Ws  is  the  weight  of  the  snow  above  the  layer  for 
which  the  density  change  is  being  computed  and  is 
expressed  in  water  equivalent  (m),  and  r?  is  the  viscos¬ 
ity  coefficient  of  snow  (rrrhr)  and  is  a  constant  for  a 
given  density,  temperature,  and  snow  type.  Kojima 
indicated  that  t?  and  ps  can  be  expressed  by 


T)  =  T?C  exp (C2ps) 


(6) 
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where  Tfc  h  the  value  of  17  when  ps  =  0,  and  C2  is  a 
constant  to  be  determined.  If  eq  5  and  6  are  com¬ 
bined  and  C1  =  eq  5  becomes 

=  C]Wiexp(-C2pi).  (7) 

Kojima  reported  C1  values  of  2.6-9.0/m  ‘hr  and  C2 
values  of  21  m3/Mg.  Equation  7  accurately  described 
the  data  of  Kojima  except  for  the  cases  of  low-density 
snow  layers,  wind-packed  snow,  and  depth-hoar  lay¬ 
ers.  For  low-density  and  wind-packed  layered  snow, 
the  density  due  to  compaction  increased  at  a  higher 
rate  than  that  predicted  using  the  C j  and  C2  values 
obtained  for  ordinary  snow.  However,  for  depth- 
hoar  layered  snow,  the  density  increase  was  at  a 
much  lower  rate. 

Based  on  his  observations  and  work  in  polar  regions, 
Mellor  (1964)  claimed  that  the  value  of  C1  varied  with 
snow  temperature  and  type.  He  provided  the  follow¬ 
ing  relation  of  r}e  (at  0°C)  to  t?0  (at  0°C): 


Yosida  (1955)  and  were  found  to  be  temperature- 
dependent.  Yosida  reported  that  the  rate  of  increase 
in  grain  diameter  at  -20°C  is  about  60%  of  that  at 
-6°C.  The  effect  of  settling  on  density  change  was 
found  to  be  important  only  in  the  early  stages  after 
snowfall.  Gunn  (1965)  indicated  that  layers  of  new 
snow  settled  at  about  1%  per  hour  immediately  after 
snowfall;  he  found  this  rate  to  be  independent  over  the 
density  range  from  0.05  to  0.15  Mg/m3.  Destructive 
metamorphism  was  found  to  be  a  slow  process  where 
ps  was  higher  than  0.25  Mg/m3 . 

There  is  no  established  mathematical  expression 
describing  density  changes  due  to  destructive  meta¬ 
morphism.  Anderson  (1976),  using  the  same  reason¬ 
ing  as  that  given  in  relating  the  density  increase  due 
to  compaction,  hypothesized  the  following  relation: 

=  c3exp[-c4(r0-ni  (10) 

tor  ps  <  pd  and 


% 

—  =  exp  - 
^0 


d  ZcT| 

R  T-T0  \ 


(8) 


1  fos 

Ps  dt 


=  C3  exp(-C4(r0-r)l  exp[-46(ps-pd)] 


where  T  and  are  the  absolute  temperatures  in  K 
that  correspond  to  6°  and  0°C,  respectively.  A  is 
the  activation  energy  104  cal/mol)  and  R  is  the 
universal  gas  constant  (~  2  cal/mol  K).  For  tempera 
tures  normally  experienced  in  areas  with  seasonal 
snow  cover,  the  value  of  A/(RTT0)  can  be  taken  as 
0.08/K  (if  T0  =  273  K  and  T  =  253  K).  Thus,  to  in¬ 
clude  the  effect  of  temperature  on  the  density  change 
due  to  compaction,  eq  7  can  be  rewritten  as 


(ID 

for  ps  >  pd.  C3  is  the  fractional  settling  rate  (hr-1 )  at 
0°C  for  ps  <  pd,  C4  is  a  settling  parameter  (K~l ),  pd 
is  the  density  below  which  the  settling  rate  for  snow 
equals  C3,  and  46  is  an  empirical  dimensional  con¬ 
stant  (m3/Mg)  necessary  to  lower  the  settling  by  a 
factor  of  100  when  p5~pd  =  0.1  Mg/m3  [that  is, 
exp(-46x0.1)  =  exp(-4.6)  ^  1/100] . 


=c,  W%  exp(-c2ps)  exp(-0.08(r0-  T)  1 


(9) 

where  C ^  is  now  the  fractional  increase  in  density 
(rrf1  hr"1 )  at  0°C  and  ps  =  0. 

Destructive  metamorphism 

Under  equilibrium  temperature  conditions,  water 
molecules  move  on  the  snow  crystals  by  the  processes 
of  sublimation  and  condensation  in  order  to  decrease 
the  surface  free  energy.  Freshly  fallen  snow  crystals 
have  a  very  high  ratio  of  surface  area  to  mass.  The 
water  molecule  migration  process  (called  destructive 
metamorphism)  changes  these  sharp-edged  crystals 
into  aggregates  of  smooth  grains  that  are  rounded, 
oblong  or  irregular.  As  a  result,  the  snow  settles  and 
increases  in  density.  These  phenomena  have  been  ob¬ 
served  and  photographed  by  Bader  et  al.  (1939)  and 


Constructive  metamorphism 

Constructive  metamorphism  is  the  process  of  vapor 
transfer  within  the  snow  cover  due  to  the  temperature 
gradient.  Vapor  is  removed  from  one  crystal  by  sub¬ 
limation  and  deposited  on  another  by  condensation. 
The  change  in  density  with  respect  to  time  is 


dPs 
3 1 


(12) 


where  m  is  the  net  sublimation  (the  net  amount  of 
vapor  that  undergoes  a  phase  change)  and  bC/bt  is 
the  rate  of  change  of  saturated  vapor  concentration 
in  the  void  space  of  snow.  When  the  effective  diffu¬ 
sion  coefficient  of  water  vapor  Dt  varies  with  snow 
depth,  the  net  sublimation  can  be  expressed  as 


m  =  De 


32C 

dz2 
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dDt  3  C 
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bt 
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where  z  is  one  of  the  coordinates  in  the  direction  of 
diffusion. 

For  water -vapor -saturated  porous  media  such  as 
snow,  the  value  of  vapor  concentration  C  is  solely  a 
function  of  temperature;  that  is,  C  =  /( 7  )  disregarding 
the  surface  energy  effect.  Substituting  ee  13  into  eq 
12  yields 


bp  37- 

— -  -  D  f'  —  +/' — +D  /" 
bt  eJ  bz2  1  bz  bz  eJ 


(14) 

where  /'  and  f"  are  the  first  and  second  partial  deriva¬ 
tives  of  vapor  concentration  with  respect  to  tempera¬ 
ture  (bC/bT  and  b2C/bT2).  If  De  is  assumed  to  be 
solely  a  function  of  temperature  for  a  given  depth 
[that  is,  bDJbz  can  be  replaced  by  (bDe/bT)(bT/bz)  \ , 
eq  14  becomes 


(15) 


high-density,  wet  snow  with  nearly  spherical  grains, 
the  rate  of  increase  in  density  should  be  similar  to  that 
for  dry  snow.  For  fresh,  low-density  snow,  it  seems 
reasonable  to  expect  that  the  presence  of  liquid  water 
will  accelerate  destructive  metamorphism,  thus  in¬ 
creasing  the  settling  rate. 


REGELATION 


The  change  in  equilibrium  melting  temperature 
dTm  due  to  a  small  change  in  hydrostatic  pressure 
is  expressed  thermodynamically  as 


(16) 


where  v  and  s  are  volume  and  entropy,  respectively,  of 
a  unit  mass  of  material,  and  subscripts  £  and  s  refer  to 
liquid  and  solid,  respectively.  Equation  16  can  be  re¬ 
written  as 


dT„ 


-  dp  -  -  A  dp 
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For  a  prolonged  thermal  gradient  near  the  bottom 
of  the  snow  cover,  constructive  metamorphism  leads 
to  the  formation  of  depth  hoar.  Kojima  (1967)  re¬ 
ported  that  well-established  depth-hoar  layers  com¬ 
pact  at  a  greatly  reduced  rate.  DeQuervain  (1973) 
also  indicated  that  layers  in  a  state  of  advanced  con¬ 
structive  metamorphism  usually  do  not  settle  unless 
there  is  structural  collapse. 

Melt  metamorphism 

Melt  metamorphism  is  the  change  in  snow  struc¬ 
ture  due  to  melt-freeze  cycles  and  the  change  in  crys¬ 
tals  due  to  the  presence  of  liquid  water.  In  general, 
melting  decreases  the  depth  of  snow  cover,  and  all  or 
a  portion  of  the  melt-water  may  be  retained  and  may 
refreeze,  causing  an  increase  in  the  ice  content  of  the 
snow  cover.  A  melt-refreeze  cycle  increases  the  den¬ 
sity  of  the  affected  portion  of  the  snow  cover  by  sev¬ 
eral  percent.  Wakahama  (1968)  reported  that  the 
grain  size  increased  at  a  faster  rate  as  the  amount  of 
water  increased,  but  the  density  (initial  ps  =  0.39 
Mg/m3)  did  not  change  unless  the  snow  was  subject¬ 
ed  to  a  load.  Colbeck  (1973)  studied  the  problem 
theoretically  and  arrived  at  essentially  the  same  con¬ 
clusions  about  grain  growth  and  density  change. 

Both  Wakahama  and  Colbeck  indicated  that  when 
water  saturation  is  high,  the  rate  of  compaction  should 
increase.  However,  these  conditions  normally  do  not 
occur  in  a  snow  cover  except  over  impermeable  ice 
layers  or  at  the  ice/snow  interface.  Therefore,  for 


where  Lf  is  latent  heat  of  fusion  and  is  defined  as 
7"m(sc-ss).  For  ice,  >  i/c;  therefore,  Tm  decreases 
with  increasing  hydrostatic  pressure. 

The  first  theoretical  value  of  A  -  0.00745°C/bar 
was  calculated  by  Thomson  (1849),  and  Thomson’s 
brother  (1850)  verified  experimentally  that  7*m  is 
lowered  by  pressure.  The  modern  accepted  experi¬ 
mental  value  of  A  (0.00738°C/bar)  was  obtained  by 
Moser  (1929).  However,  for  a  high  pressure  form  of 
ice,  >  vs;  consequently,  Tm  increases  with  increas¬ 
ing  pressure.  It  should  also  be  noted  that  under  hydro¬ 
static  tension,  ice  should  be  able  to  exist  at  tempera¬ 
tures  well  above  0°C,  but  this  has  never  been  verified 
experimentally. 

The  melting  point  is  also  changed  by  non-hydrostatic 
stress.  The  phase  equilibrium  of  a  solid  under  non¬ 
hydrostatic  stress  was  examined  theoretically  by  Gibbs 
(1 877).  He  reported  that  the  stress  component  to  be 
used  to  determine  dTm  is  the  stress  normal  to  the  inter¬ 
face.  On  the  other  hand,  Verhoogen  (1951)  indicated 
that  the  relevant  stress  is  the  mean  stress  in  the  solid. 
Kamb  (1961)  analyzed  these  two  theories  as  well  as 
the  most  recent  ones  and  concluded  that  only  that  of 
Gibbs  has  any  validity. 

The  lowering  of  7*m  of  ice  by  pressure  has  been 
used  to  explain  the  phenomenon  of  regelation,  a  term 
which  was  introduced  by  Tyndall  and  Huxley  (1857) 
to  account  for  Faraday’s  observation  that  two  pieces 
of  ice  adhere  when  they  are  brought  into  contact. 
However,  regelation  has  been  used  more  specifically 
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to  describe  the  passage  of  a  weighted  object  through 
a  block  of  ice,  the  compaction  of  wet  snow  (Colbeck 
1979),  and  the  movement  of  glaciers  over  obstacles 
on  their  bed  (Weertman  1957).  Bottomley  (1872) 
conducted  the  first  experiments  on  the  passage  of  a 
wire  through  ice  and  explained  his  observations  in 
terms  of  Thomson’s  pressure-melting  theory.  The 
pressure  of  the  wire  causes  the  ice  underneath  it  to 
melt  and  form  a  thin  layer  of  water,  which  moves 
around  the  upper  side  of  the  wire,  where  it  freezes. 

As  the  ice  continues  to  melt  and  refreeze,  the  wire 
passes  through  it.  Bottomley  concluded  that  the 
wire  has  to  be  a  sufficiently  good  thermal  conductor 
to  permit  the  latent  heat  of  fusion  needed  to  melt  the 
ice  to  be  conducted  through  the  wire  from  its  upper 
side,  where  the  water  is  refreezing  and  liberating  heat. 
Turpin  and  Warrington  (1884)  repeated  Bottomley’s 
experiment  using  wires  of  different  thermal  conduc¬ 
tivities.  They  reported  that  the  speed  of  a  wire  pass¬ 
ing  through  the  ice  increased  with  the  wire  conduc¬ 
tivity,  but  they  failed  to  develop  a  simple  relation  be¬ 
tween  the  two  parameters.  Kojima  (1954)  conducted 
a  number  of  experiments  on  regelation  using  six  wires 
with  differences  in  either  thermal  conductivity  or 
wire  diameter.*  He  concluded,  as  had  Bottomley 
(1872)  and  Turpin  and  Warrington  (1884),  that  the 
most  important  factor  for  regelation  is  the  thermal 
conductivity  of  the  wire  material  and  that  the  contri¬ 
bution  due  to  heat  conducted  along  the  wire  from 
the  surrounding  air  was  insignificant.  He  also  indi¬ 
cated  that  the  viscoelastic  property  of  ice  did  not  play 
a  role  in  the  process. 

Recently  a  series  of  theoretical  and  experimental 
papers  have  been  published  on  this  subject.  This  in¬ 
cludes  the  work  of  Telford  and  Turner  (1963),  Nye 
(1967, 1973),  Frank  (1967),  Townsend  and  Vickery 
(1967),  Nunn  and  Rowell  (1967),  Hahne  and  GrigulJ 
(1969,  1972),  and  Drake  and  Shreve  (1973).  Nye 
(1967)  analyzed  the  problem  theoretically,  taking 
into  account  the  flow  of  heat  and  water.  He  showed 
that  if  the  thermal  resistance  of  the  object  was  the 
dominant  factor,  then  the  wire  speed  depends  on  the 
volume  of  the  object  and  not  its  shape.  He  further 
reported  that  for  a  sphere  or  a  cylindrical  wire,  the 
thickness  of  the  water  film  is  always  uniform  and 
does  not  depend  on  the  speed. 

Nye  developed  the  following  equation: 

FA  A/Z.f  ypj  (18) 


♦  In  one  case  a  combination  of  high-  and  low-conductivity 
wires  was  used.  This  was  done  by  fastening  linen  cords  on 
both  ends  of  a  copper  wire  of  1 0-cm  length  and  laying  it 
over  the  ice  pillar  so  that  the  copper  wire  part  was  on  its  top. 


where  v  is  the  steady  state  velocity,  F  is  the  load  on 
the  object,  Lf  is  the  latent  heat  of  fusion,  pj  is  the  ice 
density,  X  is  the  thermal  conductivity  of  the  object  of 
volume  V,  and  A  is  a  constant  defined  in  eq  1 7.  Figure 
6  summarizes  the  results  of  Nye’s  calculations.  Metal 
wires  used  in  laboratory  experiments  usually  fall  into 
regions  1 A  or  2A;  poor  conductors  are  usually  in  1  B 
or  1C.  The  water  film  was  found  to  be  very  thin;  for 
example,  for  a  steel  wire  1  mm  in  diameter,  the  theory 
predicts  a  water  film  thickness  of  0.36  pm. 

To  test  Nye’s  theoretical  treatment,  Townsend  and 
Vickery  (1967)  conducted  experiments  with  spheres 
and  discs  that  were  about  10  mm  in  diameter  and  made 
of  different  materials.  Their  results  for  spheres  are 
shown  in  Figure  6.  The  observed  velocities  were  up  to 
three  times  lower  than  those  predicted  by  eq  18.  How¬ 
ever,  for  discs,  the  velocities  varied  from  nearly  equal 
to  the  predicted  value  to  33  times  slower.  Nunn  and 
Rowell  (1967),  using  cylindrical  specimens,  showed 
similar  discrepancies  for  metallic  wires  (by  a  factor  of 
up  to  18)  but  found  reasonable  agreement  with  the 
theory  when  testing  insulated  wires. 

Frank  (1967)  speculated  on  the  effect  of  water-film 
stability  on  the  rate  of  penetration  of  the  object.  For 
high-thermal-conductivity  objects,  the  freezing  surface 
at  the  rear  of  the  object  is  likely  to  be  unstable,  caus¬ 
ing  the  water  to  separate  from  it  in  an  extended  wake 
and  resulting  in  slower  penetration  than  the  theory 
predicts.  For  objects  with  low  thermal  conductivity, 
there  would  be  instability  in  the  forward  melting  sur¬ 
face.  However,  Frank  believed  it  would  not  produce 
any  significant  deviations  from  Nye’s  theory. 

These  studies  were  all  conducted  when  the  ice  was 
nearly  at  its  melting  point.  The  only  significant  study 
relating  the  effect  of  ice  temperature  on  the  speed  of 
wire  passage  is  reported  by  Telford  and  Turner  (1963). 
In  their  experiment,  a  steel  wire  0.45  mm  in  diameter 
was  hung  across  a  block  of  ice.  The  results  are  shown 
in  Figure  7.  The  velocity  increased  by  a  factor  of  10 
as  the  temperature  increased  from  -3.5°  to  -0.7°C. 

As  the  temperature  increased  to  -0.5°C,  the  velocity 
increased  discontinuously  by  a  factor  of  200.  This 
steep  jump  was  attributed  to  the  onset  of  pressure 
melting,  since  the  load  on  the  wire  lowered  the  melting 
point  by  about  0.5°C.  The  much  lower  velocities  at 
temperatures  well  below  -0.5°C  could  not  have  been 
due  to  pressure  melting.  Telford  and  Turner  attempt¬ 
ed  to  deduce  the  wire  motion  mechanism  at  lower 
temperatures  by  considering  the  flow  of  a  thin  New¬ 
tonian  shear  layer  of  viscous  fluid  around  the  wire. 

By  assuming  that  the  fluid  acts  as  a  liquid-like  layer 
(Fletcher  1962, 1968),  the  velocity  v  may  be  shown 
to  be  related  to  the  thickness  h  and  viscosity  rj  of 
this  layer  by 
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Figure  6.  Theoretical  and  experimental  data  on  the  passage  of  an  object  through 
ice  by  regelation,  (Adapted  from  Nye  1967,)  m  case  1,  Tr  «  l;  in  case  2,  Tr-1; 
and  in  case  3,  Tr  »  1,  where  Tr  is  the  ratio  of  the  temperature  drop  across  the  water  to  that 
across  the  object.  In  subdivision  A  the  heat  flows  mainly  through  the  object  and  the  water, 
in  B  the  heat  is  divided  between  this  path  and  that  through  the  ice,  and  in  C  the  heat  flows 
mainly  through  the  ice.  The  controlling  thermal  resistance  (ice,  object  or  water)  is  noted  in 
parentheses .  Dotted  boundaries  refer  to  spherical  objects;  all  others  refer  to  cylindrical  ob¬ 
jects.  O-  Townsend  and  Vickery's  experiments  on  spheres.  •  —  Nunn  and  Rowell's  experi¬ 
ments  on  wires. 


Figure  7.  Velocity  of  a  steel  wire  moving  through  ice  as  a 
function  of  temperature.  (Adapted  from  Telford  and  Tur¬ 
ner  1963.)  The  wire  diameter  was  0.45  mm,  the  load  was  2. 1  kg, 
and  the  thickness  of  Ice  was  10  mm.  At  the  temperature  marked 
by  the  arrow,  the  velocity  changed  discontinuousfy  by  a  factor  of 
200  due  to  pressure  melting.  The  theoretical  values  are  based  on 
Fletcher's  liquid-like  layer  theory. 


8 


V  = 


(19) 


wher$  F  is  the  force  per  unit  of  wire  length  and  a  is 
the  wire  radius.  The  variation  of  velocity  with  tem¬ 
perature  might  therefore  be  explained  by  a  change  in 
layer  thickness  with  temperature,  as  predicted  by 
Fletcher  (1962).  The  computed  values  of  v  (from 
eq  19),  with  the  thickness  and  temperature  relation¬ 
ship  given  by  Fletcher  and  the  value  of  17  appropriate 
to  supercooled  water,  are  also  shown  in  Figure  7; 
they  compare  fairly  well  with  the  experimental  data. 

Equation  19  predicts  a  linear  relation  between  v 
and  T,  but  the  experimental  results  showed  that  v  is 
proportional  to  F3,  which  is  similar  to  the  finding  by 
Glen  (1952)  for  the  creep  of  polycrystalline  ice.  How¬ 
ever,  there  are  great  differences  between  the  values  of 
activation  energy  in  Glen's  results  and  those  shown  in 
the  figure.  Therefore,  the  wire  motion  is  not  due  to  a 
simple  creep  phenomenon. 

Hahne  and  Grigull  (1972)  considered  the  regelation 
of  ice  as  a  heat  conduction  problem.  In  developing 
their  model,  they  reasoned  that  while  the  wire  is  mov¬ 
ing  through  the  ice,  the  layer  immediately  below  the 
wire  melts.  The  water  containing  the  heat  of  fusion 
is  pressed  around  the  wire  to  the  upper  side,  where  it 
refreezes  and  releases  the  heat  of  fusion,  which  is 
then  conducted  to  the  melting  zone.  For  wires  of 
thermal  conductivity  greater  than  that  of  ice,  say 
copper,  most  of  the  heat  is  conducted  through  the 
wire.  For  materials  such  as  Perlon,  however,  the  heat 
is  predominantly  conducted  through  the  ice  around 
the  wire.  The  wire  surface  temperature  should  be 
lowest  on  the  leading  side  and  highest  on  the  trailing 
side.  Hahne  and  Grigull  solved  the  heat  conduction 
equation  using  simplified  boundary  conditions  for 
the  linear  case  (when  the  heat  flows  through  the  wire) 
and  the  sinusoidal  case  (when  heat  flows  around  the 
wire).  The  penetration  rate  was  found  to  be 
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for  the  linear  thermal  boundary  and 
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volumes  of  ice  and  water,  respectively,  AP0  is  the  maxi¬ 
mum  excess  pressure,  6  is  the  water  layer  thickness,  6n 
is  the  thickness  of  the  nth  coating,  and  X,  \w,  Xj  and  Xn 
are  the  thermal  conductivity  of  the  wire  material,  the 
water,  the  ice  and  the  nth  coating,  respectively. 

To  calculate  the  penetration  rate,  the  value  of  6  is 
needed.  By  considering  the  problem  as  one  of  one¬ 
dimensional  creeping  flow,  values  of  6  can  be  approxi¬ 
mated  by  trial  and  error  from 
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for  the  linear  case  and 
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for  the  sinusoidal  case. 

Hahne  and  Grigull  used  bare  copper  and  Perlon 
wires  in  their  study  because  the  wires  differed  greatly 
in  thermal  conductivity.  The  results  can  be  repre¬ 
sented  by  a  log-log  plot  of  v  versus  a  (Fig.  8).  For 
copper  wire,  the  dependence  of  v  on  a  is  nearly  identi¬ 
cal  for  all  the  three  pressure  values  used  [v  «  (tf)'0,5)  * 
This  is  also  true  for  Perlon  wire,  except  that  the  de¬ 
pendence  of  v  on  a  is  stronger  [v «  (a)*0-1  j.  Cross 
plots  of  v  versus  pressure  p  for  any  specific  value  of 
a  used  in  the  study  revealed  the  pressure  dependence 
ofv&  (p)1-333  for  both  copper  and  Perlon  wires  (Fig. 
9).  Consequently  the  penetration  rate  can  be  repre¬ 
sented  by 

v  =  1.674  (p)1  333  (o)*0  50  (24) 

for  copper  and 

v  =  0.1 14(p)1,333  (o)“1-0  (25) 

for  Perlon  wires. 

The  effect  of  thermal  conductivity  Xmon  v  was 
evaluated  for  a  wire  radius  of  0.25  mm  and  pressures 
of  5,  7.5  and  1 0  bars  by  estimating  values  of  v  for 
Perlon  wire  from  eq  25.  It  can  be  expressed  as 


for  the  sinusoidal  case,  where  and  are  the  specific 
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Figure  8.  Penetration  velocity  as  a  function  of  wire 
radius  and  materia i  under  various  pressures.  ( From 
data  obtained  by  Hahne  and  Grigull  1972.) 


Figure  9.  Penetration  velocity  as  functions  of  pres¬ 
sure  and  wire  radius  for  two  materials  that  differ  in 
thermal  conductivities. 


v  -  0.0624(p)1-333  (Xm)0-284 .  (26) 

The  general  relationships  among  v,p}  a  and  Am  can 
not  be  established  due  to  the  lack  of  data.  However, 
the  data  available  show  that  when  a  >  0.25  mm,  the 
dependence  of  v  on  \m  is  greater;  when  a  <  0.25  mm, 
the  influence  of  Xm  on  v  is  reduced. 

The  assumption  and  verification  of  the  existence 
of  a  verv  thin  film  of  water  on  the  surface  of  ice  and 
in  any  gap  between  ice  and  any  other  material 
(Fletcher  1962,  Nakamura  1966)  led  Weyl  (1951)  to 
hypothesize  that  different  penetration  speeds  for  dif¬ 
ferent  materials  should  be  explained  on  the  basis  of 
film  thickness  rather  than  thermal  conductivity.  To 
examine  Weyl’s  reasoning,  Hahne  and  Grigull  (1972) 
compared  the  speed  of  silver-plated  copper  and  iron 
wires  with  that  of  silver  wire.  If  the  surface  material 
and  the  ice  would  dominate  the  regelation  phenome¬ 
non,  the  velocities  of  bare  silver  and  silver-plated  cop¬ 
per  wires  should  be  the  same  for  corresponding  di¬ 
ameters  and  pressures.  They  found,  however,  that 
the  same  difference  in  velocity  of  bare  copper  and 
iron  wires  is  maintained  for  the  corresponding  silver- 
plated  wires,  indicating  that  the  dominant  factor  in 
wire  movement  through  ice  is  thermal  conductivity. 
When  varnish-coated  copper  and  iron  wires  are  tested, 
the  value  of  v  and  the  difference  in  v  between  them 
are  drastically  reduced,  indicating  the  significant  ef¬ 
fect  of  wire  thermal  conductivity  on  penetration  rate. 


When  Perlon  wire  (which  has  a  thermal  conductivity 
of  the  same  order  of  magnitude  as  the  varnish-coated 
copper  and  iron  wires)  is  tested,  v  is  also  of  the  same 
order  of  magnitude.  From  these  studies,  it  is  clearly 
demonstrated  that  the  conductivity  of  the  wire  mate¬ 
rial  and  not  the  formation  of  a  thin  water  layer  be¬ 
tween  the  wire  and  ice  affects  the  wire  speed. 


Pressure ,  bars 
(•)  10 
{*)  7.5 
(o)  5.0 
(a)  8.5 


Stnusoidai  Cots 
Lineor  Cast 


0.3 

a  (mm) 


Figure  1 0.  Ratio  of  theoretical  to  experimental  pene¬ 
tration  velocities  of  Hahne  and  Grigull  (1972)  as  func 
tions  of  pressure ,  wire  radius  and  material. 
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Table  1.  Comparison  of  theoretical  and  experimental  results  on  penetration  velocity  k 


Nye /Nunn  and  Rowell _ Nye/Hahne  and  Grigull 


Materia! 

X  a  P  vtn  venr 

(W/m  K)(Wy  m)(bar)  (mm/hr)  (mm/hr) 

ytn^enr 

(mm/hr) 

yts 

(mm/hr) 

vtnlvts 

Vw 

Silver 

419 

0.25 

3.48 

120.24 

7.20 

17 

46.60 

56.41 

2.13 

6.47 

Copper 

385 

0.36 

2.41 

63.72 

3.60 

18 

24.53 

30.00 

2.12 

6.81 

0.55 

3.83 

74.88 

6.12 

12 

29.24 

35.94 

2.08 

4.78 

0.79 

5.42 

83.52 

7,20 

12 

31.81 

39.06 

2.64 

4.42 

0.79 

1.81 

27.79 

2.52 

11 

10.62 

13.04 

2.13 

4.21 

Carbon- 

50 

0,19 

8.42 

110.88 

9.00 

12 

44.48 

55.99 

1.98 

4.94 

Steel/lron 

0.23 

3.62 

39.96 

4.68 

8.6 

16.05 

20.24 

1.97 

3.42 

0.23 

7.24 

79.92 

9.72 

8.2 

32.10 

40.48 

1.97 

3.30 

0,28 

5.96 

54.72 

5.76 

9.5 

22.05 

27.84 

1.96 

3.82 

0.28 

2.98 

27.40 

2.88 

9.5 

11.02 

13.92 

1.96 

3.82 

Nylon/Perlon 

0.25 

0,21 

7.65 

5.04 

5.04 

1.0 

2.04 

2.60 

1.94 

0.41 

0.43 

3.72 

1.20 

1.73 

0.69 

0.48 

0.62 

1.94 

0.28 

0.43 

5.58 

1.80 

2.59 

0.69 

0.73 

0.93 

1.94 

0.28 

vln  =  theoretical  velocity  by  Nye  (1967). 

utr\r  ~  experimental  velocity  by  Nunn  and  Rowell  (1967). 

=  theoretical  velocity  by  Hahnc  and  Grigull  (1972)  with  linear  boundary. 
vts  =  theoretical  velocity  by  Hahne  and  Grigull  (1972)  with  sinusoidal  boundary. 


Figure  10  shows  the  comparison  between  experi¬ 
mental  results  with  bare  wires  and  Hahne  and  Grigull's 
theoretical  results  based  on  the  simplified  heat  con¬ 
duction  model.  For  copper  wire,  the  ratio  vx/vt  (where 
vx  and  ut  are  the  theoretical  and  experimental  pene¬ 
tration  velocities,  respectively)  reaches  a  value  or  ibout 
2.70  at  5.0  bars  and  a  =  0.25  mm  with  sinusoidal  boun¬ 
dary  conditions.  This  ratio  increases  as  pressure  de¬ 
creases.  However,  it  decreases  as  the  wire  radius  in¬ 
creases  for  both  the  linear  and  sinusoidal  cases.  For 
Perlon  wires,  the  ratio  for  the  sinusoidal  case  is  about 
33%  of  the  experimental  value,  and  the  effect  of  pres¬ 
sure  can  be  considered  to  be  insignificant. 

Table  1  shows  some  comparisons  of  the  theoretical 
results  of  Nye  (1967)  and  Hahne  and  Grigull  (1972) 
and  the  experimental  work  of  Nunn  and  Rowell  (1967). 
As  indicated  by  the  ratio  vtn/vts,  the  results  of  Hahne 
and  Grigull  using  a  sinusoidal  thermal  boundary  (yts) 
are  about  half  of  those  predicted  by  Nye  (vtn)  for 
both  metallic  and  insulated  wires  with  great  differ¬ 
ences  in  thermal  conductivity.  Furthermore,  the  ratio 
of  (the  penetration  rate  based  on  a  linear  thermal 
boundary)  to  yenr  (the  experimental  results  of  Nunn 
and  Rowell)  is  about  2.5  times  lower  than  the  ratio 
^tn/^enr  ^  wires.  (This  ratio  is  not  shown 

in  the  table.)  Because  of  the  different  parameters 
used  in  Hahne  and  Grigull's  experiment,  a  direct  com¬ 
parison  with  Nunn  and  Rowell's  work  can  not  be 
made.  However,  the  results  of  Hahne  and  Grigull 
were  consistently  much  higher.  Since  contamination, 
either  on  the  wire  surface  or  in  the  ice,  decreases  the 
wire  speed,  higher  speeds  are  considered  to  be  more 


reliable.  The  great  discrepancies  shown  by  vtn/venr  can 
be  attributed  to  too  high  values  of  vtn  and  too  low 
values  of  venr. 

The  most  recent  and  most  comprehensive  theoreti¬ 
cal  and  experimental  work  on  pressure  melting  and  re¬ 
gelation  of  ice  is  that  of  Drake  and  Shreve  (1973), 

They  controlled  temperature  more  precisely,  used 
much  smaller  driving  stresses,  and  watched  the  process 
and  its  after-effects  more  closely  than  did  other  investi¬ 
gators.  The  wires  were  made  of  copper,  chromel  and 
nylon  (polyamide)  with  diameters  ranging  from  0.12 
to  0.50  mm  and  Xm 's  ranging  from  0.245  to  388  W/m 
K.  They  reported  that  the  wire  speed  increases  non- 
linearly  at  all  but  the  lowest  driving  stresses;  at  a  stress 
of  about  1  bar  it  jumps  sharply  but  continuously  and 
reversibly  by  an  amount  ranging  from  six-fold  for  ny¬ 
lon  wires  to  60-fold  for  copper  wires  (Fig.  11).  Above 
this  transition  the  speeds  of  copper  wires  are  as  low  as 
one-eighth  of  those  predicted  by  Nye,  although  the 
speeds  of  nylon  and  chromel  are  about  the  same  as 
predicted.  Below  the  transition  the  speeds  of  all  wires 
are  much  less  than  predicted. 

All  wire  speeds  are  significantly  reduced  by  air  bub¬ 
bles  in  the  ice.  Drake  and  Shreve  noticed  that  the 
wires  leave  a  trace  below  the  transition  that  consists 
of  widely  scattered,  generally  tiny  bubbles  of  unfrozen 
water;  above  the  transition  the  trace  grades  from  nu¬ 
merous  bubbles  of  water  and  vapor  for  highly  conduc¬ 
tive  wires  to  a  central  flat  layer  of  water  for  poorly 
conductive  ones. 

Measurements  of  the  fractional  volume  of  water  in 
the  trace  show  that  above  the  transition,  heat  flows  to 
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Figure  1 1.  Penetration  velocities  as  functions  of  wire  mate¬ 
rial \  wire  radius  and  stress  level.  Dots  represent  the  data  of 
Drake  and  Shreve  (1973).  The  solid  lines  represent  the  the¬ 
oretical  work  of  Nye  (1967). 


the  moving  wire  from  the  surrounding  ice.  The  non¬ 
linearity  and  low  speed  below  the  transition  are  due 
to  the  accumulation  of  solutes  in  the  water  layer 
around  the  wire;  they  concentrate  toward  the  rear, 
lowering  the  freezing  temperature  and  reducing  the 
heat  flow  toward  the  melting  front.  The  transition 
occurs  when  the  temperature  at  the  rear  reaches  the 
triple  point  corresponding  to  the  existing  local  pres¬ 
sure.  With  increasing  driving  stress,  the  mean  pressure 
around  the  wire  increases  and  the  mean  temperature 
decreases,  causing  heat  to  flow  toward  the  wire  and 
causing  formation  of  the  trace,  which  carries  away 
the  volutes  For  highly  conductive  wires  (he  trace  is 
bubbly  because  of  the  Frank  instability  (1967)  of  the 
trec/mg  vurface.  which  permits  fingers  of  water  and 
vapor  to  grow  until  p  me  bed  off  by  surface  tension 
f  <>»  poorty  conductive  wires  the  non  imeanty  above 
the  transition  is  due  to  the  additional  melting 

ji  the  bon i  the  wire  and  the  ctun*  m  pressure  dis 
mhuiHin  the  •  *t  associated  with  the  forma 

imn  >>*  the  trace  f  •*  conductive  wees  the 

r*«r<  >■ near  b  9*4  i*w  me«pr\ ird  slowness  above  the 


transition  are  mainly  due  to  the  supercooling  required 
for  a  finite  rate  of  freezing,  which,  like  the  presence  of 
dissolved  solutes,  lowers  the  freezing  temperature  at 
the  rear  of  the  wire. 

The  work  of  Tozuka  et  al.  (1979)  is  the  most  re¬ 
cent  addition  to  the  literature  on  this  subject.  In  their 
experiment,  copper  and  nylon  wires  of  various  radii 
were  used  and  experiments  using  both  commercial  ice 
and  ice  doped  with  hydrogen  fluoride  were  conducted 
in  a  room  maintained  at  1°C;  the  pressure  ranged  from 
2  to  15  bars.  They  observed  a  sharp  jump  in  copper 
wire  speed  (about  1 .5  times)  at  pressures  around  5  bars. 
(This  is  in  sharp  contrast  to  the  60-fold  jump  in  wire 
speed  at  pressures  around  1  bar  reported  by  Drake  and 
Shreve  (1973),]  In  the  case  of  nylon  wires,  however, 
the  transition  pressure  has  never  been  observed.  (Drake 
and  Shreve  observed  a  6-fold  increase  in  speed  for  Per- 
lon  wires.) 

The  Findings  of  Tozuka  et  al.  concerning  the  effect 
of  wire  radius  on  speed  were  essentially  consistent  with 
the  data  obtained  by  Hahne  and  Grigull  (1972).  For 
low  thermal  conductivity  wires  (nylon  and  Perl  on),  the 
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wire  speeds  were  found  to  be  inversely  proportional 
to  the  wire  radius  (1  la);  for  copper  wires,  the  speeds 
were  proportional  to  (I/0)0*3  [versus (l/o)0'5  derived 
from  Hahne  and  Grigull’s  data], 

Tozuka  et  at.  also  studied  the  effect  of  thermal 
conductivity  on  wire  speed.  At  pressures  below  the 
transition  point  (such  as  2.4  bars),  measured  speeds 
per  unit  pressure  (v/p)  were  much  smaller  than  those 
expected  from  the  theory  (Nye  1967)  and  had  a  larger 
scatter  regardless  of  thermal  conductivity.  For  pres¬ 
sures  greater  than  the  transition  (such  as  8  bars), 
speeds  of  medium  thermal  conductivity  wires  such 
as  chromel  and  constantan  were  about  the  same  as 
predicted  from  the  theory.  But  in  the  case  of  copper 
wire,  with  its  much  higher  thermal  conductivity  value, 
the  theoretical  wire  speed  was  three  times  larger  than 
the  measured  one.  On  the  other  hand,  the  measured 
speed  for  nylon  wires,  which  are  poor  conductors, 
was  greater  than  the  theoretical  one  by  about  40%. 

The  effects  of  impurity  on  wire  speed  were  also 
reported  by  Tozuka  et  al.  They  indicated  that  the 
abrupt  change  in  wire  speed  at  the  transition  was 
mainly  caused  by  impurities  in  the  ice.  They  demon¬ 
strated  this  by  measuring  the  temperature  difference 
between  the  top  and  bottom  sides  of  the  wire. 

Regardless  of  the  uncertainties,  however,  the  ex¬ 
perimental  and  theoretical  results  as  a  whole  clearly 
demonstrate  that  the  long-accepted  explanation  of 
the  motion  of  wires  through  ice  as  a  process  of  pres¬ 
sure  melting  and  regelation  is  basically  correct.  The 
large  qualitative  and  quantitative  discrepancies  be¬ 
tween  the  simple  theory  and  the  experimental  obser¬ 
vations  are  largely  due  to  the  neglect  of  the  dissolved 
solutes  in  the  water  layer  around  the  wire,  of  the  for¬ 
mation  of  a  trace  behind  the  wire,  whose  geometry 
is  governed  by  the  Frank  instability,  and  of  the  super¬ 
cooling  required  to  freeze  water  at  a  finite  rate. 

THERMAL  PROPERTIES  OF  SNOW 
AND  FRESH  WATER  ICE 

Heat  capacity  of  snow  and  ice 

The  heat  capacity  is  defined  as  the  heat  required 
for  a  unit  of  mass  to  rise  one  unit  of  temperature  at 
constant  pressure  cp  or  at  constant  volume  cv.  Mathe¬ 
matically  it  is  defined  as 


where  H  and  U are  enthalpy  and  internal  energy  per 
unit  mass,  respectively.  Since  the  heat  needed  to 
warm  up  the  air  and  vapor  in  the  interstices  is  very 
small,  the  heat  capacity  of  dry  snow  and  ice  are  essen¬ 


tially  equal.  The  value  of  cv  can  be  computed  from 
by 


where  V  is  the  volume  of  ice,  7C  is  the  coefficient  of 
volumetric  expansion,  and  to  is  the  compressibility. 

The  value  of  cv  is  about  3%  less  than  cp  at  the  melting 
point,  and  the  difference  (cp  -cv)  decreases  with  de¬ 
creasing  temperature. 

The  values  of  the  heat  capacities  of  ice  made  in  the 
early  twentieth  century  were  reviewed  and  summarized 
by  Dickinson  and  Osborne  (1915)  and  Dorsey  (1940). 
Giauque  and  Stout  (1936),  in  an  extensive  study,  re¬ 
ported  cp  values  of  hexagonal  ice  over  temperatures 
ranging  from  1 5  to  273  K.  Flubacher  et  al.  (1 960) 
obtained  cp  values  for  hexagonal  ice  in  the  range  from 
2  to  27  K.  Sugisaki  et  al.  (1968)  also  made  extensive 
measurements  of  cp  for  amorphous,  cubic  and  hexa¬ 
gonal  ice  between  20  and  250  K.  Most  recently,  Ash¬ 
worth  (1972),  using  a  new  technique,  determined 
values  of  cp  between  50  and  270  K.  Figure  12  shows 
all  the  experimental  data  on  cp  for  hexagonal  and 
cubic  ice.  The  data  were  divided  into  three  or  four 
temperature  ranges  [i.e.,  15-50  K,  50-95  K  (or  15-95 
K),  95-150  Kand  >  150  K)  and  subjected  to  linear 
regression  analysis.  For  each  range  cp  is  expressed  as 
cp  =  A+BT.  The  constants  A  and  B  and  the  correla¬ 
tion  coefficients  0  are  listed  in  Table  2. 

The  data  from  all  the  investigations  are  consistent, 
as  the  high  values  of  0  indicate.  The  effect  of  axis- 
orientation  is  within  the  experimental  error.  The 
linear  representation  of  the  data  can  be  used  to  deter¬ 
mine  the  value  of  cp  at  any  specific  temperature.  For 
example,  cp  at  the  normal  melting  point  is  37.7  J/mol 
K,  about  one-half  of  the  heat  capacity  of  liquid  water. 
The  value  of  cp  decreases  with  decreasing  T  and,  as 
expected,  approaches  zero  at  0  K. 

Anderson  (1976)  gives  an  expression  for  as 

1.6738  +  0.1327  7*  (29) 


Table  2.  Constants  A  and  B  in  cp  =4  +  flTand0. 


Temp. 

range 

(K) 

A 

B 

<t> 

15-50 

-2.2171 

0.2094 

0.9920 

50-95 

-0.0035 

0.1597 

0.9892 

15-95 

-0.8994 

0.1710 

0.9931 

95-150 

2.2841 

0.1  3S0 

0.9817 

>150 

2.7442 

0.1282 

0.9737 
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Figure  1 2.  Specific  heat  of  ice.  Points  represent  the  data  of  Giauque  and 
Stout  (1936),  Flubacher  et  al.  (I960),  Sugisaki  et  a!.  (1968)  and  Ashworth 
(19 72).  Samples  consisted  of  cubic  and  hexagonal  ice. 


where  c\  is  the  specific  heat  of  ice.  At  193  K,  eq  29 
shows  Cj  to  be  0.2  J/mol  K  less  than  the  value  com¬ 
puted  from  the  highest  temperature  range  in  Table  2. 
At  273  K,  however,  it  shows  Cj  to  be  0.055  J/mol  K 
higher. 

Latent  heat 

The  change  in  enthalpy  dH  associated  with  a  change 
in  phase  of  a  material  at  constant  pressure  p  is  given  by 

dH  =  dU+pdV  (30) 

where  dU  and  dV  are  the  changes  in  internal  energy 
and  volume  per  unit  mass  of  material,  respectively. 

The  latent  heat  of  fusion  of  ice  Lf  is  defined  as  the 
change  in  enthalpy  when  a  unit  mass  of  ice  is  con¬ 
verted  isothermally  and  reversibly  into  liquid  water. 
Because  the  term  pdV  is  rather  small  in  comparison 


to  other  terms,  Lf  is  nearly  equal  to  the  change  in  in¬ 
ternal  energy.  Measurements  of  Lf  made  prior  to  1925 
have  been  reviewed  by  Smith  (1925)  and  those  up  to 
1940  by  Dorsey  (1940).  Many  early  measurements  are 
in  e*ror  because  the  effect  of  impurities  in  ice  was  not 
taken  into  account  properly.  The  most  reliable  value 
for  Lf  at  0°C  and  standard  atmospheric  pressure  is  33.5 
kj/kg,  reported  by  Rossini  et  al.  (1952).  The  value  of 
Lf  decreases  linearly  with  decreasing  temperature  down 
to  about  -10°C  and  then  decreases  at  a  lower  rate  as 
the  temperature  decreases  further  (Fig.  13). 

The  latent  heat  of  sublimation  of  ice  L%  represents 
the  enthalpy  change  when  a  unit  mass  of  ice  is  trans¬ 
formed  isothermally  and  reversibly  into  water  vapor 
at  standard  atmospheric  pressure  and  at  the  ice-1  iqu  id  - 
vapor  triple  point  (273.16  K).  Rossini  et  al.  reported 
a  value  of  2838  kj/kg.  It  is  interesting  that  Lf  repre¬ 
sents  only  about  1 2%  of  Ls;  therefore,  only  about  1 2% 
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Figure  13.  Latent  heat  of  ice  as  a  function  of 
temperature. 


Figure  14.  Thermal  conductivity  of  ice  as  a  function 
of  temperature.  Points  represent  the  data  of  Jakob 
and  Erk  (1929),  Powell  (1958),  Ratcliffe  (1962),  Dean 
and  Timmerhaus  (1963),  Wolfe  and  Thieme  (1964), 
Dillard  and  Timmerhaus  (1966)  and  Ashworth  (1972). 


of  the  hydrogen  bonds  must  break  when  ice  melts. 
According  to  the  data  available,  the  value  of  L%  re¬ 
mains  practically  constant  (2838  k J /kg)  between 
213.16  and  272.16  K. 

Thermal  conductivity  of  ice 

Thermal  conductivity  X  of  a  solid  material  is  de¬ 
fined  as  the  proportionality  constant  in  the  one- 
dimensional  form  of  Fourier’s  law  of  heat  conduction; 
that  is, 


(31) 


where  q  is  the  heat  flux  in  the  direction  normal  to  the 
temperature  gradient  ddjdz.  Thermal  conductivity 
usually  depends  on  temperature  Q  and  to  a  lesser  ex¬ 
tent  on  ice  crystallographic  orientation.  Early  work 
on  thermal  conductivity  of  ice  X}  was  reviewed  by 
Fowell  (1958).  At  temperatures  near  273  K,  the  value 
of  Xj  is  about  2.2  W/m  K,  about  four  times  as 
large  as  the  Xw  of  water  at  273  K.  The  increase  in  Xj 
with  decreasing  temperature  reported  by  Lee  (1905) 
was  found  to  be  much  less  rapid  than  that  obtained  by 
Jakob  and  Erk  (1929).  Landauerand  Plumb  (1956) 
reported  no  significant  differences  in  the  thermal  con¬ 
ductivity  coefficients  of  laboratory-grown  single  crys¬ 
tals,  glacial  single  crystals  and  polycrystalline  commer¬ 
cial  ice,  although  the  Xj  along  the  C-ax/s  of  the  single 
crystals  appeared  to  be  about  5%  greater  than  that  nor¬ 
mal  to  the  C-axis. 

Figure  14  summarizes  the  recent  measurements  of 
X,  of  polycrystalline  ice  by  Ratcliffe  (1962,  taken  from 
his  best-fitted  lines),  Dillard  and  Timmerhaus  (1966), 
Ashworth  (1972),  Dean  and  Timmerhaus  (1963),  and 
Wolfe  and  Thieme  (1964).  The  relatively  old  data  of 
Jakob  and  Erk  (1929)  and  Powell  (1958)  are  also 
shown.  (For  clarity  of  presentation,  only  one  notation 
is  used.)  It  seems  as  if  each  investigator  provided  data 
under  a  different  set  of  techniques  and  conditions 
(such  as  sample  purity,  sample  preparation,  and  meth¬ 
ods  for  obtaining  X).  Since  there  was  a  lack  of  data 
between  150  and  195  K,  the  data  were  divided  at  that 
gap.  Regression  analyses  were  made  based  on  the 
expression 

Xj  =  j  exp(hF).  (32) 

Table  3  shows  the  values  of  a,  b  and  <f>.  The  values 
for  the  whole  temperature  range  had  the  highest  correla¬ 
tion  coefficient.  Therefore,  for  practical  purposes,  it  is 
suggested  that  one  use 

Xj  =  9.828  exp(-0.0057  T)  (33) 


to  compute  the  values  of  Xj.  The  scattering  of  data  is 


Table  3.  Values  of  a,  b  and  0  from  the 
regression  analysis. 

Temperature 

(K) _ a _ b _ 0 _ 

<150  12.285  -0.0076  0.7970 

>195  6.727  -0.0041  0.5962 

all  T  9.828  -0.0057  0.9313 


mostly  due  to  differences  in  sample  purity  and  prepara¬ 
tion,  adjustment  of  heat  loss  to  the  surroundings,  and 
reproducibility  of  experimental  data. 

Thermal  conductivity  of  snow 

The  heat  transfer  processes  in  snow  are  much  more 
complicated  than  in  ice.  In  snow,  heat  is  transferred 
by  conduction  through  the  interconnected  ice  grains, 
by  conduction,  convection  and  radiation  across  the 
air  space,  and  by  the  movement  of  vapor  by  sublima¬ 
tion  and  condensation.  In  the  determination  of  ther¬ 
mal  conductivity,  a  temperature  gradient  is  imposed, 
which  subsequently  establishes  a  vapor  gradient  and 
thus  causes  vapor  diffusion.  Therefore,  the  thermal 
conductivity  of  snow  includes  vapor  diffusion.  Ther¬ 
mal  conductivity  is  expressed  as  effective  thermal 
conductivity  Xse  (Yen  1962)  to  account  for  all  the 
processes  occurring  in  the  snowpack.  However,  be¬ 
cause  of  the  low  temperatures,  the  effect  of  radiation 
transfer  is  usually  not  significant. 

Numerous  investigators  have  reported  the  values  of 
XS€,  and  without  exception,  they  empirically  correlated 
their  results  with  snow  density  (Fig.  15)  as  the  sole  par¬ 
ameter.  The  temperature  ranges  under  which  the  data 
were  obtained  were  usually  not  defined.  The  scatter¬ 
ing  of  data  is  due  to  the  snow  conditions  (such  as  aging 
and  grain  size  distribution)  and  to  the  effects  of  vapor 
diffusion.  Another  factor  may  be  the  method  used  to 
determine  X^.  The  transient  method  is  more  accurate 
than  the  steady-state  technique,  because  in  the  transi¬ 
ent  method  the  heat  loss  to  the  surroundings  during 
the  experiment  need  not  be  assessed.  Regardless  of 
the  difference  in  results  from  each  investigator,  all 
data  can  be  reasonably  represented  by 

X„  =  2.22362  (p$)1885  (34) 

with  a  0  of  0.861 4.  It  gives  a  reasonable  value  of  X{ 
when  extrapolated  but  a  lower  value  of  Xa  when 
Pa- 

Pitman  and  Zuckerman  (1967)  appear  to  be  the 
only  investigators  who  systematically  considered  the 
effect  of  temperature  on  Xje.  They  used  vapor-grown 
ice  crystals  and  conducted  Xse  measurements  at  -5°, 
-27°,  and  -88°  C  and  at  densities  ranging  from  0.1  to 


0.6  Mg/m3.  A  semi-log  plot  (Fig.  16)  of  their  data 
(taken  from  their  graph)  shows  three  parallel  straight 
lines  (with  the  exception  of  ps  =  0.1  Mg/m3  at  -88°C). 
Since  the  thermal  conductivity  of  air  at  0 °C  is  about 
0.024  W/m  K,  the  straight  line  extrapolation  to  ps 
=  0.1  is  probably  a  better  representation  of  their  data 
at  -88°C.  To  determine  the  effect  of  temperature  on 
Xsc  from  these  data,  a  plot  of  Xse  (taken  from  Figure 
16  at  any  specific  ps)  versus  0  was  made  and  Xse  was 
found  to  be  proportional  to  0;  that  is,  Xse  a  exp 
(0.00880).  Finally,  the  data  of  Pitman  and  Zucker¬ 
man  can  be  represented  [except  for  the  data  for  the 
lowest  density  (ps  =  0.1 )  ]  by 

Xse  =  0.0688  exp(O.OO880 +  4.6682  ps)  (35) 

as  shown  in  Figure  1 7.  Also  shown  in  the  figure  are  a 
few  points  taken  from  the  work  of  Yen  (1965)  and 
Schwerdtfeger  (1963a).  For  ps  -*pj  =  0.917  and 


p%  (Mg/m3) 

Figure  15.  Effective  thermal  conductivity 
of  snow  as  a  function  of  density.  Poin  ts 
are  derived  from  the  equations  of  Abel’s  (1893), 
X^  =  2.8451 (o$)\  Jansson  (1901),  X^  =  0.0209 
+  0. 7950(ps)7  +  2.5104(ps)* ;  Devaux  (1933), 

\se  =  0.0293 + 2. 9288(p$)7 ;  Kondratevo  (1945), 
X*  *  3.5564(ps)i ;  Bracht  (1949),  X^  *  2.0502 
(ps)*;  Sufakvelidee  (1959),  X^  *  0.51 05  p; 

Yen  (1962),  X^  *  3.221 7(ps). 


16 


0  0.2  0.4  0.6 

f>%  (Mg/m3) 


Figure  /  6.  Effective  thermal  conductivity  of  snow  Figure  1 7.  Effect  of  temperature  and  snow  density  on 

as  func  tions  of  sno  w  density  an  d  temperature.  e  ffecti ve  thermal  conducti vity  of  sno  w. 


6  =  0°C,  the  value  of  A$e  A-,  =  4.98  W/m  K,  which  is 
about  2.3  times  the  accepted  ice  value  of  2.2  W/m  K 
at  0°C.  For  ps  -  pd  and  0  =  0°C,  Ase  -+  Aa  =  0.0688 
W/m  K,  which  is  about  2.8  times  the  value  for  air  at 
0°C  (Xa  =  0.0247  W/m  K). 

Schwerdtfeger  (1963a),  using  Maxwell’s  work  on 
the  electrical  conductivity  of  heterogeneous  media, 
derived  theoretical  expressions  for  Xse.  For  dense 
snow,  X5e  is  equivalent  to  Xja  and  is  given  as 


where  5  is  a  constant  related  to  porosity  by  r?  =  1/(1  +s)3. 
When  tj  is  given,  the  value  of  5  can  be  calculated.  When 
ps  =  0.4  and  77  =  0.536, 5  is  found  to  be  0.21 2.  From 
eq  38,  Asa  =  0.7023  W/m  K  for  snow  at  0°C.  For  very 
low  density  snow,  a  model  is  developed  consisting  of  a 
suspension  of  small  spherical  particles  in  air  not  in  con¬ 
tact  with  each  other;  the  model  is  identical  to  eq  36 
except  for  the  interchange  of  A|  and  Aa ,  and  the  substi¬ 
tution  of  Asa  for  Aia. 


2Xj+Xa-2i}(Xj-Xa) 

Aii"  2x,+xa-n(xrxa)  ' 


where  Xja  is  the  thermal  conductivity  of  dense  snow 
or  ice  containing  air  bubbles,  and  rj  is  the  porosity  and 
is  related  to  pt  and  ps  by  r?  =  1  -(p/p,).  Because  Aa 
«  X;  and  n  =  1  -(pj/ps),  eq  36  becomes 

‘-■i&v  13,1 

For  low  density  snow  (down  to  0.1  S  Mg/m3) 


.  (2*5)5  ■ 

(1+J)2 


Effective  thermal  diffusivity 

For  a  homogeneous,  isotropic  medium  with  con¬ 
stant  physical  properties,  the  effective  thermal  diffusiv¬ 
ity  ae  is  the  coefficient  of  the  general  heat  conduction 
in  Cartesian  coordinates;  that  is, 


where  at  is  defined  as  Ae/psc j  and  psCj  is  the  volumetric 
specific  heat  of  the  medium. 

Sulakvelidze  (1959)  formulated  a  heat  transfer  equa¬ 
tion  for  porous  media  containing  saturated  vapor,  water 
or  ice  at  temperatures  close  to  those  of  phase  transition. 
This  was  done  by  including  an  evaporation-condensation 


17 


term  in  Fourier’s  heat  conduction  equation.  For  the 
case  of  one-dimensional  heat  transfer,  the  heat  con¬ 
duction  equation  becomes 


can  be  considered  to  be  constant.  In  one-dimensional 
steady  diffusion  processes,  Fick’s  law  gives  the  mass 
flux/77  as 


bT 

bt 


£l 

bz2 


(40) 


-  ndC 

m  -  -  D  — 
dz 


(46) 


where  q  is  the  intensity  of  sublimation  of  condensa¬ 
tion  (Mg/m3 )  and  can  be  expressed  in  terms  of  vapor 
concentration  C  by 


d-De 


92C 

dz2 


dC 

dr 


(41) 


where  Dt  is  the  effective  water  vapor  diffusivity 
(m2/s).  Substitution  of  a  prescribed  function  f(T) 
for  C  in  eq  41  yields 


<7=Oe 


if  _2£ 

dz2  dt 


(42) 


Equations  40  and  42  can  be  combined  to  give 
[1+OeLs/(aeCiPs)/'l  d2y 


bT  t-s^aecjP$)j  i 

d7=a«  1  HLJc,pAf’  f 


(Wcipjr  (dT\2 


1  +(LJcxp%)  f  ('3// 


(43) 


where  /'  and  /"  are  first  and  second  derivatives  of  / 
with  respect  to  T.  Equation  43  describes  the  general 
processes  of  heat  transfer  taking  place  in  a  moist  por¬ 
ous  medium  without  including  the  effects  of  solar 
radiation  or  the  changes  in  the  amount  of  liquid 
water  within  the  snow  medium.  Since  the  values  of 
Cjps  are  much  greater  than  those  of  LJ\ eq  43  can 
be  reduced  to 


dT/  D'LS\&T  ,  WV-* 
dt  \°e  cipi)dz2  *  cipi  \dzi 


(44) 


When  there  is  no  vapor  diffusion  (£s  =  0),  eq  4  re¬ 
duces  to  the  well-known  Fourier  equation 


d  T  d2T 

*  "  dz2  ‘ 


(45) 


Once  the  function  /  is  chosen,  the  temperature  distri¬ 
bution  as  a  function  of  time  can  be  computed. 


Heat  transfer  by  water  vapor 
diffusion  in  snow 

The  diffusion  coefficient  D  for  a  binary  system  is 
a  function  of  temperature,  pressure  and  composition, 
For  low  pressure  gas  mixtures  or  dilute  solutions,  D 


where  dC/dz  is  the  concentration  gradient  in  the  di¬ 
rection  of  diffusion. 

Vosjda  (1950)  was  the  first  to  study  the  effective 
diffusion  coefficient  De  for  water  vapor  diffusing 
through  snow.  He  reported  Dt  values  ranging  from 
0.7  to  1  .Ox  1 CT4  m2/s,  which  is  about  four  or  five 
times  larger  than  the  diffusion  coefficient  Dq  for 
water  vapor  diffusing  through  air.  When  p$  is  in  the 
range  of  0,08  to  0.5  Mg/m3,  Dt  values  remain  more  or 
less  constant.  By  maintaining  the  lower  boundary  tem¬ 
perature  higher  than  the  upper  one,  he  found  that  the 
effect  of  natural  convection  is  insignificant.  To  explain 
the  fact  that  De  >  DQl  Vosida  pointed  out  that  ice 
grains  do  not  act  as  mere  obstacles  to  the  diffusion  of 
water  vapor  as  sand  grains  do;  they  produce  water 
vapor  themselves,  thereby  facilitating  the  mass  trans¬ 
fer  processes  of  sublimation  and  condensation.  The 
heat  transfer  due  to  molecular  vapor  diffusion  can  be 
written  in  a  form  similar  to  Fourier’s  law  of  heat  con¬ 
duction: 


(47) 


where  qv  is  the  heat  flux  due  to  vapor  diffusion  and 
0j  is  the  ratio  of  water  vapor  density  to  temperature. 
The  product  &jDeLs  can  be  considered  as  Xv,  the  ther¬ 
mal  conductivity  due  to  molecular  diffusion.  If  Ls 
=  2828.38  k J/kg,  0T  =  0.39x  1 0"3  kg/m3  K  and  Oe 
=  0.85x  1( T4  m2/s,  then  \v  =  0.094  W/m  K.  There¬ 
fore,  the  molecular  vapor  diffusion  plays  a  significant 
role  in  the  process  of  heat  transfer  in  low  density  snow. 


Heat  and  vapor  transfer  with 
forced  convection 

The  effects  of  air  flow  through  snow  on  Ae  and  De 
have  been  studied  experimentally  and  theoretically 
and  were  reported  in  a  series  of  papers  by  Yen  (1962, 
1963,  1965).  For  one-dimensional  and  steady-state 
cases,  the  governing  equations  for  heat  and  vapor  trans¬ 
fer  can  be  summarized  as  follows.  For  heat  transfer, 


dT  GMwLtdP  d2 T 

Gc>j^—z^r-^+\~=0  <48) 


Mn  dz 


dz 4 


and  for  vapor  transfer, 
d2P , 


tlM  *vr"  *h\dJ!±  -n  MQ1 

dz2  W  D'tdPJdz  dt  /  dz  °  (  9* 
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where  ps  is  the  saturation  vapor  pressure  of  snow,  ca 
is  the  heat  capacity  of  air,  /?v  is  the  gas  law  constant 
for  water  vapor,  G  is  the  air  mass  flow  rate,  Mw  and 
M  are  the  molecular  weights  of  water  and  dry  air,  re* 
spectively,  7r  is  the  total  pressure  of  the  system, 
is  the  mean  temperature  of  the  system,  A^  and  D'e  arc 
the  effective  thermal  conductivity  and  diffusivity 
under  forced  convective  flow,  respectively,  and  0  is 
defined  as  0  =  GR^TmMw/MTrD'e.  To  calculate  A^, 
only  the  steady-state  temperature  distribution  is  need¬ 
ed;  to  obtain  ,  the  density  distribution  before  and 
after  the  experiment  and  the  steady-state  temperature 
distribution  are  required. 

For  unconsolidated  snow  with  a  density  between 
0.376  and  0.472  Mg/m3 ,  the  value  of  A^  can  be  ex¬ 
pressed  by  Ag  =  Ae  +  24.27G  (Yen  1962),  where  Ae  =  0.586 
W/m  K  and  G  varies  from  10x10“3  to  40x1 (T3  kg/m 2  s. 
The  value  of  0.586  for  Ae  is  consistent  with  the  data  re¬ 
ported  by  Kondrateva  (1945).  The  value  of  D\  can  be 
represented  by  D'e  =  3.01 6x  10“ 3  (G  +  0.456x  10“3),//2 
when  G  =  0  and  D'e  Oe  =  0.65x1 0-4  m2/s.  As  indi¬ 
cated  in  the  previous  section,  Yosida  (1950),  using  a 
completely  different  approach,  reported  a  value  of 
0.85x1 0'4  m2/s  for  De.  The  differences  may  be  par¬ 
tially  due  to  the  lower  temperature  in  Yen's  experiment. 

For  naturally  compacted  snow  (Yen  1965),  the  fol¬ 
lowing  correlations  for  A^  and  D*e  are  valid  for  ps  rang¬ 
ing  from  0.50  to  0.59  Mg/m3  and  G  varying  from  5x 
HT3  to  32x1  O'3  kg/m2  s. 

X;  =  3.22(ps)2  +  25.1G  (50) 

and 

D'e  =  0.65  x  1 0’4  +  0.05 1 3(ps) 3-20  (G)061 5 

(51) 

where  and  D'e  are  expressed  in  W/m  K  and  m2/s, 
respectively.  When  no  air  is  flowing,  \'z  reduces  to 
Ae,  which  is  3.22(ps)2,  and  D’t  becomes  £e,  which  is 
0.65x1 0"4  m2/s.  Both  values  are  in  agreement  with 
experimental  results  (Yen  1962,  1963). 

These  results  show  that  air  flow  has  a  considerable 
effect  on  the  values  of  Aj,  and  D\  of  both  unconsoli¬ 
dated  and  naturally  compacted  snow.  The  increase  in 
thermal  conductivity  and  vapor  diffusivity  due  to  air 
flow  is  responsible  for  the  small  variations  in  snow 
density  and  temperature  gradient  near  the  surface 
layers  of  a  snow  cover. 

THERMAL  PROPERTIES  OF  SEA  ICE 

Compared  with  fresh-water  ice,  whose  physical 
properties  are  well  known,  sea  ice  is  a  relatively  com¬ 
plex  substance  whose  transformation  to  a  completely 


solid  mixture  of  pure  ice  and  solid  salts  is  attained 
only  at  low  temperatures  so  extreme  that  they  are 
rarely  encountered  in  nature.  The  physical  properties 
of  sea  ice  thus  depend  strongly  on  salinity,  temperature 
and  age.  Many  of  these  properties  are  still  not  fully  un¬ 
derstood,  particularly  those  important  for  the  under¬ 
standing  of  natural  ice  covers.  In  fact  it  appears  to  be 
safer  to  rely  on  theoretical  values  for  thermal  conduc¬ 
tivity  and  specific  heat,  because  precise  measurement 
of  these  properties  has  always  posed  considerable  dif¬ 
ficulty.  In  the  development  of  a  suitable  sea-ice  model 
for  calculating  these  quantities,  an  interesting  progres¬ 
sive  complexity  appears:  To  calculate  specific  heat, 
knowledge  of  the  composition  alone  is  sufficient.  Cal¬ 
culations  of  density  also  require  consideration  of  air 
bubble  content.  Finally,  to  calculate  the  thermal  con¬ 
ductivity,  information  on  the  spatial  distribution  of  a)) 
components  is  necessary. 

Specific  heat  of  sea  ice 

When  sea  water  is  cooled  to  its  freezing  point,  pure 
ice  crystals  form.  As  freezing  progresses,  pockets  of 
brine  are  cut  off,  so  that  the  sea  ice  is  composed  of 
pure  ice,  brine,  solid  salt  crystals  and  air  bubbles 
(which  have  a  negligible  effect  on  the  specific  heat). 

The  equilibrium  salt  concentration  of  the  brine  trapped 
in  the  ice  depends  on  its  temperature.  When  the  sea  ice 
temperature  increases,  the  brine  is  diluted  to  a  new 
equilibrium  concentration  by  the  melting  of  pure  ice 
at  the  ice/brine  interface.  The  specific  heat  of  sea  ice 
is  the  total  of  the  heat  required  to  raise  the  sea  ice  con¬ 
stituents  (i.e.  the  pure  ice  and  the  brine)  a  unit  of  tem¬ 
perature  and  the  heat  associated  with  the  phase  change. 
The  latter  process  leads  to  the  fact  that  sea  ice  has  an 
abnormally  large  specific  heat. 

Salinity  is  usually  defined  by  oceanographers  as  the 
number  of  grams  of  dissolved  solids  in  one  kilogram  of 
solution.  This  also  applies  to  sea  ice,  and  the  mass  of 
salt  in  grams  per  kilogram  of  sea  ice  is  usually  quoted 
in  parts  per  thousand  (%o)«  However,  for  convenience 
in  deriving  phase  compositions  of  sea  ice,  the  units 
“grams  of  salt  per  gram  of  sea  ice"  are  often  used. 

Based  on  the  work  of  Malm^ren  (1927),  it  is  believed 
that  salinity  is  usually  a  good  parameter  in  analyzing 
the  thermal  behavior  of  sea  ice. 

Specific  heat  of  ice  between  the 
freezing  point  and  -8.2° C 

Within  this  temperature  range,  there  is  no  signifi¬ 
cant  substitution  of  ions  in  the  ice  lattice;  thus,  the 
brine  in  the  interior  cells  has  the  same  salt  composition 
as  the  sea  water  from  which  the  ice  was  formed.  No 
precipitation  of  salt  occurs  in  this  temperature  range 
(Fig.  18),  and  the  specific  heat  depends  essentially  on 
the  relative  extent  of  the  phase  change  and  the  specific 
heat  of  pure  ice  and  the  salt  solution.  The  thermal 
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Figure  18 .  Freezing  point  of  brine  as  a  function  of 
fractional  salt  content  (A  fter  Schwerc/tfeger  1 963b.) 


Figure  19.  Specific  heat  of  sea  ice  as  func¬ 
tions  of  temperature  and  salinity.  (Modified 
from  Ono  1966.) 


capacity  of  salt  in  4°/oo  sea  ice  is  on  the  order  of 
3.3472  J/kg  K  and  can  be  considered  to  be  negligible. 
The  effects  of  heat  of  crystallization  (or  dilution) 
can  also  be  neglected  based  on  data  quoted  by  Lange 
and  Forker  (1952).  Consequently  the  specific  heat 
of  sea  ice  csi  can  be  written  as 

csi  =  +  f  (Cw"ci)+ci  <52) 

where  a  is  the  coefficient  of  the  relation  between 
fractional  salt  content  of  the  brine  s  (in  grams  of  salt 
per  gram  of  water)  and  the  temperature  0,  o  is  the 
salinity  of  sea  ice  (grams  of  salt  per  gram  of  sea  ice), 
and  Cj  and  cw  are  the  specific  heat  of  pure  ice  and 
water,  respectively.  The  relation  between  s  and  0  is 
linear  (Fig.  18),  so  s  =  ad  ( a  is  the  slope  of  s  versus  0, 
a  negative  quantity). 

Replacings  in  eq  52,  Pounder  (1965)  and  Schwerdt- 
feger  (1963b)  both  gave  the  expression 

C*i  =  -^2^+^(C--C')  +  C*  (53) 

in  which  the  term  oCj  has  been  omitted  for  the  case 
of  natural  sea  ice.  Ono  (1966)  used  a  somewhat  dif¬ 
ferent  approach  and  developed  the  following  expression: 


csi  =  0.505  +0.001 80  +431 1 .5  -O.8o+O.2Oa0. 

(54) 


The  values  of  csj  are  shown  in  Figure  19.  Since  the 
contributions  of  the  4th  and  5th  terms  are  negligible, 
eq  54  reduces  to 


csi  =  0.505+0.00180+431 1.5  —  .  (55) 

6 2 

Specific  heat  of  sea  ice  between 
-8.?  and -23°  C 

The  continuous  deposition  of  Na2S04  •  10H2O 
takes  place  only  within  this  temperature  range  (Fig. 
18).  If  the  rate  of  deposition  with  temperature  is  lin¬ 
ear,  then  the  quantity  of  precipitation  at  any  tempera¬ 
ture  can  be  calculated  by  extrapolating  from  the  initial 
section  of  the  phase  graph.  Therefore,  in  1  gram  of  sea 
ice,  there  are  w  grams  of  water,  ws  grams  of  dissolved 
salt,  wp  grams  of  precipitated  salt  with  which  0wp 
grams  of  water  are  combined  in  the  crystals,  where  0 
is  the  mass  ratio  of  10H2O  to  Na2SO*  and  equals 
1 .27.  Thus,  in  general,  the  mass  of  pure  water  in  the 
unfrozen  solution  is 


(56) 
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where  p  is  grams  of  precipitate  (not  including  its  water 
of  crystallization). 

The  mass  of  pure  ice  is 

/»i  =  1-o-^(1+0P).  (57) 

The  change  in  mass  of  unfrozen  water  w  as  fractional 
salt  contents  and  associated  precipitate p  change  is 

dw  =  _2i*±M'  (58) 

(5+P)2 

(The  salinity  o  is  a  constant;  only  the  redistribution  of 
salt  from  liquid  to  solid  or  vice  versa  is  of  concern.) 
Figure  18  shows  that 

s+p  =  a0  and  p~  a  (0  +  8.2).  (59) 

Substitutings+p  =  ad  into  eq  58  gives  dw  =  -(o/aO2) 
dd.  Therefore,  the  heat  absorbed  by  a  unit  mass  of 
sea  ice  for  a  temperature  increase  dd  is 

dq  =  Lfdw+  |l  -a- (1  +/3p;]Ci<yfl 

+  5^C^+[^(1+0))C^  (60) 

where  ch  is  the  specific  heat  of  the  precipitated 
hydrate. 

The  expression  for  specific  heat  of  sea  ice  is  ob¬ 
tained  by  evaluating  dq/dd  of  eq  60,  substituting  eq 
59  in  the  resulting  expression  and  rearranging,  and  is 
expressed  as 

C*i=-^^  +  ^(C-'Ci)  +  (1‘a)C* 

frr0*flc„1  (61) 

by  Schwerdtfeger  and  Pounder  (1962).  Since  aja 
<  1 ,  (0  +  8.2/0)  <  1 ,  the  terms  of  0c  j  and  (1  +  0)ch  are 
of  the  order  of  1 ,  and  o  ^  1 0”3,  the  last  term  and 
0C|  can  be  neglected.  Equation  61  reduces  to  the 
same  form  shown  in  eq  52.  Values  of  csj  for  salinities 
ranging  from  0  to  0,01  gram  of  salt  per  gram  of  sea 
ice  and  temperatures  ranging  from  -2°  to  -22°C  are 
shown  in  Figure  1 9. 

Heat  of  fusion  of  sea  ice  when 
0°  >  0  >  -8.2° C 

As  sea  ice  has  no  fixed  temperature  for  phase 
transition,  the  heat  of  fusion  L  J(  of  sea  ice  is  the 


amount  of  heat  required  to  melt  1  gram  of  sea  ice  of 
salinity  a  and  temperature  0.  The  vaiue  of  LSI  is 
found  by  integrating  eq  4  from  0  to  0m#  where  0m  is 
the  temperature  when  the  melting  is  completed  and  is 
computed  from  the  relations:  mx  -  1  -a- (a/s).  When 
mK  =  0  (that  is,  when  melting  is  complete)  and  when 
5  =  a0m  is  substituted,  0m  is  0m  *  o/a,  after  taking 
into  account  that  o  «  1 .  If  higher  order  terms  of  o 
are  neglected, 

LS)  =  79.68  -  0.5050  -27.3a+431 1 .5  j 

+  O.8o0-O.OO902.  (62) 

The  values  of  LSi  from  eq  62  are  shown  in  Figure  20. 

For  practical  purposes,  the  fifth  and  sixth  terms  in 
the  right  side  of  eq  62  can  be  neglected,  and  it 
becomes 

Z.si  =  79.68-O.5O50-27.3o+4311.5|.  (63) 

Density  and  thermal  conductivity  of  sea  ice 

The  thermal  conductivity  of  sea  ice  is  strongly  de¬ 
pendent  on  composition,  that  is,  the  density,  salinity 
and  temperature.  Schwerdtfeger  (1963b)  showed  that 
thermal  conductivity  is  mainly  determined  by  salinity 
at  high  temperature  and  by  density  at  low  temperat-res. 
As  with  specific  heat,  it  is  simplest  to  consider  ice  in  the 
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Figure  20.  Latent  heat  of  fusion  of  sea  ice 
as  functions  of  temperature  and  salinity. 
(Modified  from  Ono  7966.) 
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temperature  range  between  0°  and  -*8.2°C,  below 
which  the  precipitation  of  Na2S04  *10H2O  compli¬ 
cates  the  analysis. 

Composition  and  air  bubble  content 
of  sea  ice  above  -8.2°  C 

With  a  and  s  defined  as  before,  the  mass  of  unfro¬ 
zen  brine  per  unit  mass  of  sea  ice  ismb  =  a/s+o.  The 
volume  of  brine  is  (o+o/s)/pw(1+s)  =  a/spw,  and  the 
volume  of  ice  is  (1  -o-o/s)/pi#  where  pw  and  p{  are 
the  densities  of  pure  water  and  ice,  respectively.  The 
volumes  of  brine  and  ice  per  unit  volume  of  sea  ice 
are 


Vb=p»spw 


and 


(64) 


where  psj  is  the  density  of  sea  ice.  The  volume  of  air 
per  unit  volume  of  sea  ice  is 


6  ( °c ) 


=  1  -Psi 


1  -a-a/s\ 

P\  /' 


(66) 


lfpw  =0.999  Mg/in3,pj  =  0.917  Mg/m3, 5  =  aO,  and 
a  =  - 1 .82x10"J/°C,  eq  66  reduces  to 


,  /  1-Q 

*Mo.917 


+ 


4.98  a\ 

e  r 


(67) 


Figure  21  shows  the  air  bubble  content  by  volume  of 
sea  ice  as  functions  of  temperature,  salinity  and 
density. 

Rewriting  eq  66,  replacings  =  ad ,  and  neglecting 
the  terms  containing  aaB  or  as  gives  the  sea  ice  density 
as 


P$i~  Ot0pw-ff(pw-Pj)  ' 


(68) 


Since  the  second  term  in  the  denominator  is  rather 
small,  rewriting  and  substituting  the  values  of  pw ,  pj 
and  a  changes  eq  68  to 


Psi  =  (1  ■  j1  +  tt0pw  (Pw  -P|)]pi 

=  (1-^)(l-^)°.917.  (69) 


Figure  27.  Fractional  air  content  (by  volume)  as  a 
function  of  temperature  for  sea  ice  of  different  salini¬ 
ties  and  densities. 


if  the  air  content  is  constant.  VMiile  the  density  is  tem¬ 
perature-dependent,  especially  near  the  freezing  point, 
the  dependence  diminishes  considerably  at  lower 
temperatures. 

Thermal  conductivity  model  for  sea  ice 

Langleben  (1960)  and  others  have  shown  sea  ice 
to  consist  of  pure  ice  enclosing  vertical  cylinders  of 
brine  that  are  approximately  elliptical  in  cross  section 
and  whose  lengths,  especially  at  higher  temperatures, 
are  large  compared  with  their  average  diameters.  Ander¬ 
son  (1958),  on  the  other  hand,  preferred  to  calculate 
the  thermal  conductivity  as  if  the  brine  pockets  were 
spherical.  This  assumption  may  have  greater  validity 
at  lower  temperatures,  when  the  small  amount  of 
brine  plays  a  less  significant  role  in  determining  the 
thermal  conductivity.  Based  on  Maxwell's  principle 
that  sea  ice  consists  of  uniformly  and  randomly  dis¬ 
tributed  spherical  air  bubbles,  the  thermal  conduc¬ 
tivity  of  bubble  ice  Xb{  (Schwerdtfeger  1963b)  is 

2Xj+Xa  -  2l/.(Xi-X.) 

Xbi=  2\,+X1+^(XrX1)  Xl'  (70) 


This  equation  shows  that  p$j  increases  with  salinity 


For  a  temperature  range  of  0°  to  -20°C,  and  taking  the 


Figure  22.  Density  and  thermal  conductivity  of  Figure  23.  Effective  thermal  conductivity  of  sea  ice  as 

bubbly  ice  as  a  function  of  fractional  air  content.  a  function  of  temperature  for  various  salinities  and 

densities. 


thermal  conductivity  of  ice  X;  to  be  2.09  W/m  K  and 
the  thermal  conductivity  of  air  Xa  to  be  2.5 1x1 0"2 
W/m  K,  cq  70  is  used  to  calculate  Xbj  and  is  shown, 
along  with  the  density  of  bubbly  fresh-water  icc,  in 
Figure  22. 

If  the  sea  ice  consists  of  parallel  configurations  of 
bubbly  pure  ice  and  enclosed  brine,  then  the  thermal 
conductivity  of  sea  ice  Xsi  is 

op^ 

K  =  (71) 


Thermal  diffusivity  of  sea  ice 

Thermal  diffusivity  is  the  most  directly  observable 
thermal  property,  as  it  is  directly  related  to  the  rate  of 
temperature  change  in  a  medium.  It  has  been  seen  that 
the  values  of  csj ,  psj  and  X$j  are  monotonic  functions 
of  the  temperature;  that  is,  as  temperature  increases, 
cs |  and  psj  also  increase  but  Xsi  decreases.  The  values 
of  asj  =  ^Si/Psicsi  are  sh°wn  'n  Figure  24  as  functions 
of  temperature  and  salinity.  Note  that  the  values  of 
<*sj  are  not  greatly  affected  by  density  changes  as  they 
are  for  freshwater  ice. 


where  Xb  is  the  thermal  conductivity  of  the  brine, 
which  is  strongly  dependent  on  the  concentration 
and  to  a  lesser  extent  on  the  temperature.  If  the 
change  in  thermal  conductivity  and  fractional  salt  con¬ 
tent  is  linear  at  constant  temperature  for  s  <  0.1 5,  then 
on  the  basis  of  thermal  conductivity  measurements  for 
NaCI  and  Na2S04  solutions  given  by  Lange  and  Forker 
(1952),  the  value  of  Xb  can  be  approximated  by 


Method  of  determining  thermal  diffusivity 

In  general,  the  equation  for  thermal  diffusion  is 
written  as 


30  =  ^si  30  3^0 

CsiP*i  3/  3 z  dz  Asi  dz2 

where  t  is  time  and  /  is  the  depth.  Let 


(73) 


Xb  =  0.41 84(1. 25 +0.0300 +0.0001 402).  (72) 

Using  Xb  from  this  expression,  X,  =  2.09  W/m  K.  The 
values  of  Xsj  calculated  from  eq  71  are  shown  in  Fig¬ 
ure  23. 

All  the  curves  (that  is,  for  all  densities  and  salinities) 
exhibit  asymptotic  behavior.  At  lower  temperatures, 
the  thermal  conductivity  of  low-salinity  sea  ice  is  equal 
to  the  value  for  fresh-water  ice. 


j_  8xsi/  be/dz  \ 

xsi  dz  '32,0/3/2' 


Equation  73  may  then  be  rewritten 


bd 

bt 


0+/K, 


(74) 


(75) 


Letting  0Z ,  be  the  temperature  at  depth  z  and  time  t, 


o-  (gm  solt/gm  sea  ice) 
0  x  i<53 


lOxiO*7) 


m)*o 


a.  Temperature  profile  at  times 
i) (O,.,  =  0J»  '■>(0/.s=0,.Jand 

h(o,  =  W 


Figure  24.  Thermal  diffusivity  of  sea  ice  as 
a  function  of  temperature  for  various  salini¬ 
ties.  (The  effect  of  density  on  a$j  for  a  giv¬ 
en  salinity  is  insignificant.) 


'il  1  b.  Graphic  determination 

;  ; _ j _  of  times  l/f  i  ?  and  t  ?  from 

!i  %  *s  1  actual  temperature  record. 

Figure  25.  Determining  thermal  diffusivity  from  temperature 
profiles.  (At  t2,  the  minimum  temperature  is  close  to  depth  i\ 
therefore ,  it  may  be  assumed  that  f  -  Oat  approximately  t  =  l2 . 


expressing  0Z+St,  0z_st,and  0zt+T  in  the  Taylor 
series,  and  neglecting  the  higher  order  terms  yields 


^z, t+r  *  5  (^z+s,t  +^z-s,t) 

d0  s2  d20 
~Tdf6bz2- 

From  eq  75  and  76, 


®z,t+r  =  g  (®z+s,t  +  ^®z,«  +  ®J-S,t)-  (78) 

The  quantity  on  the  right  side  of  eq  78  can  be  com¬ 
puted  from  the  experimental  data  on  the  temperature 
changes  at  depth  z-5,  z  and  z+s,  and  thus  the  value 
of  0Z  t+y,  from  which  the  value  of  r  can  be  determined 
graphically.  Using  the  values  of  5  and  r,  the  quantity 


s2/6t  can  be  computed;  it  equals  ((1  +/)arsijzt.  If  f 
=  0,  the  value  of  asj  can  be  determined  immediately. 

By  eq  74,  /  =  0  when  36 /dz  =  0,  that  is,  when  the  tem¬ 
perature  is  at  the  minimum  or  maximum  value  on  the 
vertical  profile.  Figure  25a  shows  temperature  profiles 
at  times  (0Z_S  =  0Z),  /2(0Z_S  =  0Z+S)  and  t3(0t  =  0Z+S). 
These  times  are  easily  determined  from  the  temperatures 
as  shown  in  Figure  25b.  At  each  of  times  ,  f2  and  f3, 
the  values  of  */6  (^z_s^ t  +  40z#t  +  0Z+S>t)  can  be  com¬ 
puted,  and  the  corresponding  values  of  r( ,  r2  and  r3 
can  be  determined  graphically.  With  the  values  of  r’s 
and  a  given  distance  increment  (which  may  be  varied), 
the  corresponding  values  of  [(1  +/)<*sjJZfl  can  be  ob¬ 
tained  using  eq  77.  However,  at  time  r2  the  minimum 
temperature  is  closer  to  depth  z;  therefore,  it  may  be 
assumed  that  /  =  0  at  approximately  r2,  and  the  value 
of  [(1  ’f/)otSjJz  t  reduces  to  asi  directly. 


SUMMARY 

This  report  dealt  comprehensively  with  the  thermal 
properties  of  snow,  ice  and  sea  ice.  Studies  on  ice 
density,  thermal  expansion  and  compressibility  were 
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reported.  The  available  data  were  graphed  to  make 
them  easier  to  use  and  to  demonstrate  the  large  dis¬ 
crepancies  in  the  compressibility  data.  The  mechan¬ 
isms  and  processes  associated  with  snow  density 
changes  due  to  compaction  and  destructive,  construc¬ 
tive  and  melt  metamorphism  were  discussed  in  detail. 
Equations  describing  these  processes  were  derived. 

Various  theories  of  rcgelaiion  were  reviewed.  Data 
on  wire  speed  as  functions  of  wire  material  (for  ex¬ 
ample,  high  and  low  thermal  conductivity  or  metallic 
and  non-mctallic  wires),  wire  radius,  driving  stress, 
temperature  and  contamination  were  analyzed.  Tem¬ 
perature  measurements  near  the  moving  wire  pro¬ 
duced  some  insight  into  the  processes  of  melting  and 
refree/ing.  There  is  still  no  theory  that  incorporates 
ail  the  factors  that  affect  the  rate  of  wire  penetration 
into  the  ice. 

The  determination  of  heat  capacity  and  latent  heat 
of  fusion  of  fresh-water  ice  was  well  documented. 
Slight  variations  in  the  reported  heat  capacity  values 
were  noted. 

For  values  of  thermal  conductivity  of  ice  and 
snow,  however,  a  much  greater  scatter  was  observed. 
For  ice,  the  large  variation  was  believed  to  be  due  to 
the  condition  of  the  sample  (such  as  the  level  of  con¬ 
tamination  and  the  air  bubble  content).  Reproducible 
data  for  snow  were  much  more  difficult  to  obtain. 
Because  snow  undergoes  constant  metamorphism  as 
soon  as  it  falls  on  the  ground,  its  thermal  conductiv¬ 
ity  is  evidently  not  solely  a  function  of  its  density. 
However,  thermal  conductivity  of  snow  was  usually 
reported  in  terms  of  density  alone  without  any  indi¬ 
cation  of  age  or  temperature.  (The  temperature  effect 
is  considered  to  be  small.)  However,  the  age  of  the 
snow  sample  and  the  temperature  level  under  which 
it  had  been  stored  apparently  contributed  to  the  dis¬ 
persion  and  uncertainty  in  the  data.  A  brief  review 
and  discussion  of  the  effect  of  air  ventilation  on  ef¬ 
fective  thermal  conductivity  and  vapor  diff usivity 
were  also  presented. 

Determining  the  thermal  properties  >f  sea  ice  is 
considerably  more  complicated  than  determining 
those  of  fresh-water  ice.  All  the  thermal  properties, 
such  as  specific  heat,  latent  heat  and  thermal  conduc¬ 
tivity,  arc  functions  of  both  temperature  and  salinity 
(in  the  temperature  range  0°C  >  0  >  -8.2°C).  In 
general,  the  specific  heat  is  much  higher  than  that  of 
the  fresh-water  ice,  especially  when  the  temperature 
0  -*■  0°C.  This  unusual  phenomenon  is  caused  by  the 
fact  that  when  the  sea  icc  temperature  rises,  melting 
occurs,  adding  the  latent  heat  of  fusion  to  the  sensible 
heat.  The  latent  heat  of  fusion  of  sea  icc  is,  however, 
generally  lower  than  that  of  the  fresh  water  ice. 

In  the  temperature  range  0°C  ^  0  >  - 8.2°C ,  ex¬ 
pressions  of  the  volumetric  traction  of  air  and  sea  icc 


density  in  terms  of  salinity  and  temperature  were  de¬ 
rived.  For  temperatures  lower  than  -8.2°C  the  analy¬ 
sis  becomes  much  more  complex  because  of  salt  pre¬ 
cipitation.  By  considering  sea  ice  to  be  bubbly  pure 
ice  with  brine  pockets,  and  by  assuming  that  the  com¬ 
ponents  are  arranged  in  parallel  (like  the  connection  in 
resistances),  thermal  conductivity  of  sea  ice  was  ex¬ 
pressed  in  terms  of  bubbly  pure  ice,  brine,  salinity, 
densities  of  sea  ice  and  water,  and  temperature.  Due 
to  its  complex  composition  and  its  transient  behavior 
(the  migration  of  brine  or  the  aging  effect),  sea  ice  has 
not  been  studied  extensively;  only  a  handful  of  field 
studies  have  been  reported.  In  the  conclusion  of  thi* 
report,  an  innovative  method  for  determining  the  ther¬ 
mal  d if f usivity  of  sea  ice  was  briefly  described. 
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